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NOMENCLATURE 

A„  coefficient  of  Boltzmann  equation  defined  by 

-^  Eq.  (A-18) 

B,  coefficient  of  Boltzmann  equation  defined  bv 
Zyn  Eq.  (A-18) 

C.  coefficient  of  Boltzmann  equation  defined  by 
Xn  Eq.  (A-19) 

D  coefficient  of  Boltzmann  equation  defined  by 

-^  Eq.  (A-19) 

D  diffusion  coefficient  (cm) 

fg  angular  flux  (neutrons/sec- cm2) 

H  height  of  rod-moderator  region  (cm) 

I0  (r)  ordinary  Bessel  function  of  imaginary  argument 

J0  (r)  ordinary  Bessel  function 

J  neutron  current  (neutrons/sec-cra2) 

P„  i£l)  '  associated  Legendre  polynomials 

r  distance  from  the  axis  of  the  cylinder 

Rl  radius  of  absorber  (inches) 

R2  outer  radius  of  moderator  region  (inches) 

s  displacement  in  the  direction  of  £L 

s  expansion  coefficient  defined  by  Eq.  (A-ll) 

S*  source  of  thermal  neutrons  (neutrons/sec-cm  ) 

S  expansion  coefficient  defined  by  Eq.  (A-13) 

JLm 

U©(<xr)  a  particular  solution  of  Bessel »s  equation 
(see  Appendix  C) 

o{  separations  equation  constant  defined  by 
Eq.  (26) 

A  separations  equation  constant  defined  by 

'  Eq.  (16) 


X         separations  equation  constant  defined  by 
Eq.  (26) 

V  gradient  operator 

V  Laplacian  operator 

A         separations  equation  constant  defined  by 
Eq.  (16) 

Arr        transport  mean  free  path  (cm) 

^2        ratio  of  macroscopic  absorption  cross  section 
to  diffusion  coefficient  (cm*2) 

_Q         unit  vector  in  the  direction  of  the  neutron 
motion 

<P  angle  between  the  projection  of  jQ_  on  the  x-y 

plane  and  the  direction  of  r 

<fi  thermal  neutron  flux  (neutrons/sec-cm2) 

q-  standard  deviation 

Z!  total  macroscopic  cross  section  (era"*-1-) 

£  macroscopic  absorption  cross  section  (cm""  ) 

Zs  macroscopic  scattering  cross  section  (cm  ) 

£\  macroscopic  transport  cross  section  (cm""l) 

&  angle  between  jQ.   and  the  z  axis 

Q  angle  between  the  unit  vectors  Q.     and  _Q 

Subscripts 

JL  primary  indices  in  the  spherical  harmonic 

expansion 

yn.  secondary  indices  in  the  spherical  harmonic 

expansion 

1  absorber  region 

2  moderator  region 


INTRODUCTION 

Within  the  active  core  volume,  the  various  materials 
present  compete  with  each  other  in  the  absorption  of  the  neu- 
trons present.   In  reactors  that  are  assumed  to  be  homoge- 
neous, the  ratios  of  the  cross  sections  may  be  used  to  well 
describe  the  events  that  occur.  However,  in  a  reactor  that 
is  considered  hetrogeneous,  which  is  the  usual  case,  the  neu- 
tron flux,  and  hence  the  rate  of  absorptions  and  collisions, 
may  vary  greatly  between  the  materials  present. 

Consider  the  case  where  an  absorbing  rod  is  separate 
from  the  moderating  material.  If  an  experimental  determi- 
nation is  made  of  the  neutron  flux  in  the  moderator  and  the 
rod,  a  depression  of  the  flux  within  the  rod  will  be  evident. 
This  depression,  though  expected,  is  greater  than  that  pre- 
dicted by  a  diffusion  theory  calculation  using  the  value  of 
the  flux  at  the  surface  of  the  rod  as  the  average  flux  inside 
the  absorber.  This  phenomenon  is  called  self- shielding  or 
self- absorption,  and  it  is  due  to  the  absorption  of  a  consid- 
erable fraction  of  the  incoming  neutrons  by  the  nuclei  near 
the  surface  of  the  absorbing  rod. 

Some  work  has  been  done  which  treats  self-shielding  and 
capture  of  neutrons  by  cylindrical  shapes  surrounded  by  a 
scattering  medium  (I4.,  8,  9).  Bartels  (1)  surveyed  all  the 
earlier  literature  (pre-1950)  on  self-shielding  for  various 
absorber  geometries.  He  presented  this  information  in 
graphical  form  along  with  the  analytical  expressions  used 


to  construct  the  graphs.   In  each  case,  the  problem  of  self- 
shielding  was  approached  from  the  viewpoint  of  a  self-absorp- 
tion factor  and  an  effective  cross  section.   Dworlc,  et.  al. 
(i;)  investigated  self-shielding  for  an  infinitely  long,  hollow 
cylinder.   This  group  rigorously  derived  equations  expressing 
the  self-shielding  in  a  hollow  cylinder,  and  presented  their 
results  in  graphical  form.  However,  they  too  used  a  self- 
shielding  factor  and  an  effective  cross  section  to  account  for 
the  increased  flux  depression.  Kushneriuk  (8)  investigated 
the  problem  of  neutron  capture  by  long  cylinders.  He  arrived 
at  an  expression  for  the  self-shielding  correction  in  terms  of 
the  average  flux  in  the  rod  and  the  unperturbed  flux.   So,  in 
essence,  an  effective  cross  section  was  determined. 

Thus,  it  appears  that  the  conventional  method  of  treating 
self-shielding  is  by  determining  an  effective  cross  section 
for  the  absorber.  Then,  the  proper  neutron  absorption  is 
obtained  by  multiplying  the  average  neutron  flux  inside  the 
absorber  by  the  effective  cross  section. 

A  different  approach  to  the  problem  of  self- shielding 
is  one  that  relates  the  measured  value  of  the  cross  section 
and  the  physical  radius  of  the  rod  to  an  effective  rod  radius. 
That  is,  given  a  rod  of  known  radius  and  cross  section,  an 
effective  rod  radius  is  determined  so  that  the  neutron  absorp- 
tion as  calculated  using  the  effective  rod  radius,  the  neutron 
flux  at  the  surface  of  the  rod,  and  the  cross  section,  would 
match  the  experimentally  determined  flux  depression. 


Thus,  to  fully  utilize  the  concept  of  an  effective  rod 
radius,  one  would  like  to  have  available  graphs  that  relate 
the  effective  rod  radius  to  some  dimensionless  parameter  that 
is  peculiar  to  the  system  under  consideration.   Since  the 
product  of  the  macroscopic  absorption  cross  section  and  the 
rod  radius  is  dimensionless,  graphs  involving  this  quantity 
as  well  as  moderator  properties  would  reduce  the  labor 
involved  in  the  calculation  of  the  thermal  neutron  flux  level 
and  the  power  distribution  throughout  the  active  core  volume. 


THEORY 

A  rigorous  solution  of  the  Boltzmann  transport  equation 
using  spherical  harmonics  results  in  an  expression  for  the 
flux  that  is  an  infinite  series  involving  the  associated 
spherical  harmonics,  i.e.,  the  P*   (jQ-  )  or  the  associated 
Legendre  polynomials.   By  considering  only  a  finite  number 
of  terms  in  the  series,  an  approximation  to  the  value  of  the 
flux  can  be  made.   The  accuracy  of  the  approximation  depends 
on  the  number  of  terms  of  the  infinite  series  that  are  speci- 
fied.  These  approximations  are  symbolically  represented  as 
P  approximations,  where  the  value  of  JL   denotes  the  order  of 
the  approximation;  e.g.,  JL   =  1  is  a  first  order  approximation, 
A    =  2  is  a  second  order  approximation,  and  so  on.   The  sub- 
script ~?n  ranges  from  -JL   to -*-^  .  A  more  detailed  discussion 
of  the  spherical  harmonics  solution  of  the  Boltzmann  trans- 
port equation  in  cylindrical  geometry  is  presented  in  Appendix 
A,  while  a  general  treatment  of  this  subject  can  be  found  in 
references  (3,  5,  7,  11,  16). 

The  P^  approximation  results  when  the  infinite  series 

representation  of  the  flux  is  terminated  after  two  terms, 

J.    -   0  and  J,  -   1.   The  three  equations  that  are  obtained  lead 

to  a  solution  for  the  flux  that  is  identical  with  the  result 

of  ordinary  diffusion  theory  (3,  7,  12,  16)  except  for  the 

definition  of  D,  the  diffusion  coefficient.   For  the  ?± 

approximation,  the  diffusion  coefficient  is  defined  as 

l/(3E,  ),  while  for  diffusion  theory,  the  diffusion  coeffi- 
tr 


cient  is  defined  as  1/(3 Zs  ).  Keeping  in  mind  this  differ- 
ence in  notation,  the  mathematical  methods  associated  with 
diffusion  theory  will  be  applied  to  the  problem  under  consid- 
eration. 

The  geometry  of  this  problem  is  shown  in  Figure  1.  The 
center  region,  region  I,  represents  an  absorber.  The  z  axis 
is  co-linear  with  the  axis  of  region  I.  Region  II,  which  is 
concentric  to  region  I,  represents  the  moderator.  The  source 
is  located  in  the  x-y  plane,  and  it  can  assume  any  configura- 
tion, i.e.,  point  source,  plane  source,  discrete  sources,  and 
so  forth. 
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Figure  1.  Coordinate  System 
As  neutrons  from  the  source  diffuse  in  the  positive  z 
direction,  a  depression  of  the  neutron  flux  in  region  I,  the 
absorber  region,  is  apparent  when  the  flux  is  plotted  as  a 
function  of  r,  the  distance  from  the  z  axis,  on  any  given  z 
plane.  The  depression  is  caused  by  the  excess  absorption  of 
neutrons  in  region  I  over  that  of  region  II.  However,  the 


magnitude  of  the  flux  depression  due  to  these  absorptions  is 
greater  than  would  be  predicted  by  diffusion  theory. 

The  increase  in  the  flux  depression  is  the  result  of 
self-absorption  within  the  absorbing  media.  That  is,  the 
incident  neutron  flux  will  be  decreased  by  absorptions  at  the 
surface  of  the  media  which  decreases  the  number  of  neutrons 
available  for  absorption  as  the  neutrons  diffuse  inward  through 
the  absorber.  To  obtain  the  true  absorption  rate,  the  flux 
could  be  determined  as  a  function  of  r,  multiplied  by  the 
macroscopic  absorption  cross  section,  and  the  product  inte- 
grated over  the  volume  of  the  absorber.  This  method  is  cumber- 
some and  difficult  to  use  (1).  A  simpler  procedure  would  be 
to  devise  an  effective  cross  section  to  multiply  the  flux  at 
the  surface  of  the  absorber  to  yield  the  correct  neutron  absorp- 
tion (k) •     A  third  method  of  obtaining  the  proper  neutron 
absorption,  and  the  method  of  this  work,  is  by  increasing  the 
radius  of  a  fictitious  absorbing  media  while  allowing  the 
macroscopic  absorption  cross  section  of  this  media  to  remain 
its  measured  value. 

Thus,  for  an  absorber  with  a  given  macroscopic  absorption 
cross  section,  an  effective  absorber  radius  can  be  calculated 
by  varying  the  radius  of  the  fictitious  absorber  until  the 
function  that  best  describes  the  flux  within  this  absorber 
fits  the  measured  values  of  the  flux  with  the  minimum  least 
squared  error.  Keeping  the  above  in  mind,  the  mathematical 
theory  of  this  problem  will  now  be  presented. 


Using  the  subscript  "1"  for  the  absorber  and  the 
subscript  "2"  for  the  moderator,  the  equations  to  be 
satisfied  are 

and 


(1) 


Vz<k    -ycA~  O  (2) 


'Z'Z 


where 


and 


xf«*sA*;  **  =  WM 


0^  =  1/(31  );        D,=  V<3Z   ). 


(3) 


(i+) 


In  order  to  obtain  a  unique  solution  to  Eqs.  (1) 
and  (2),  boundary  conditions  must  be  applied.  For  this 
type  of  problem,  the  boundary  conditions  are  well  known 
(3»  5*  6,  7,   11 j  12)  and  straight  forward.  The  flux,  4>   , 
must  be  finite  and  non-negative  throughout  the  region  of 
diffusion.  This  condition  is  clear  since  the  flux  can 
never  be  infinite  or  negative.  However,  it  can  be  zero. 
At  interfaces  between  media  with  different  diffusion 
properties,  the  net  neutron  current  densities  and  the 
neutron  fluxes  are  equal.  Finally,  the  neutron  flux 
must  go  to  zero  at  the  extrapolated  boundary.  The  extrap- 
olated boundary  is  defined  as  the  physical  boundary  plus 
0.71A  ,  where,  in  general,  A+r  is  small  compared  to  the 
dimension  of  the  physical  boundary. 


8 


The  following  five  equations  will  express  mathemat- 
ically the  above  boundary  conditions: 

1.  4j(0,y)  ¥  oo  (5) 

2.  <tHRlty)~0  (6) 

3.  <^M)  =  ^(^H)  =0  (7) 
k.          4><{R\,y)   =  ^(R^y)  (8) 

The  boundary  condition  at  the  x-y  plane  depends  upon  the 
source  configuration.  For  this  problem  it  will  be  assumed 
that  the  source  is  a  plane  isotropic  source  emitting  Sf 
neutrons  per  second  per  centimeter  squared.  Thus,  the 
following  two  boundary  conditions  are  obtained: 

b<h 


6.    S'/Z    =5=  -D, 


7.  S'/Z     =  S  =  -D,4^ 


-o 


=  o 


(10) 


(11) 


The  Laplacian  operator,  V   $   for  cylindrical 
coordinates  with  no  G    dependence  is  written  as 

Vi=fe  +7TF  +  -$".  (12) 

Using  this,  Eq.  (1)  for  region  I  becomes 

?&  +  ±A±  +£±l  -  it**  =o.  (13) 


If  the  variables  are  separable,  a  solution  of  the  form 
is  attempted.   Utilizing  Eq.  (li|)  in  Eq.  (13)  produces 


U'^R'+ltf-tf.O 


(15) 


R  x   rR 

Because  the  first  two  terms  are  functions  of  r  alone, 
they  must  be  equal  to  a  constant,  A    .  By  similar 
reasoning,  the  third  term  is  equal  to  a  constant,  A    . 
Thus  Eq.  (15)  becomes 

/32  +  A2-XN  0  (16) 

and  this  equation  will  be  called  the  separations  con- 
stants equation. 

Now,  examining  the  individual  terms  of  Eq.  (15)» 
it  is  seen  that 

&  R"  +  rfc  *  - /S*  (17) 

and 

j  Tl   -  A*.  (18) 

A  solution  of  Eq.  (17)  is 

R(r)  =  AXa(/?r;  +  BK0(y3rJ  (19) 

while  a  solution  for  Eq.  (16)  is 

1(y)  -  C  SmK  Ay  +  E  Cosh  Ay.  (20) 
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Thus,  the  complete  solution  before  the  application  of 
the  boundary  conditions  is 

^f»;^)=[AIor/6r)+BK.(/3r)][CSmh  Aj,  +  ECosb  A^].   (21) 

In  order  that  Eq.  (5)  be  satisfied 

B=0.  (22) 

Changing  the  argument  of  Eq.  (20)  from  Ky   to  A(H~3), 
which  is  valid  for  a  slab  geometry  coordinate,  Eq,  (7) 
is  satisfied  if 

E=  0.  (23) 

Thus, 

<^irty)  =  AI0(/Sr)SmK  \{\\-y)  .  (24) 

Applying  the  same  procedures  to  region  II,  the 
solution  to  Eq.  (2)  may  be  written  as 

<jt(r;^*[F0o(^r)  +  GYo(^r)][MSmhy^+  NCosfctf^]     (25) 

where  the  separations  constants  equation  for  this  region 
is 

-ex*   +  V7-  -  «l   =   0,  (26) 

The  choice  of  sign  on  the  ck1  term  is  determined  by  the 
type  of  solution  desired  in  the  outer  region.  Applica- 
tion of  Eqs.  (6)  and  (7)  to  Eq.  (25)  produces 
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<kbtf=  8rj0(o<r)-  TY0(«r)]Sinh  WH-jJ  (27) 

where 

T=  J0(<*R2)/Y0(*<R2)  •  (28) 

Now,  an  auxiliary  function  UQ(  <X   r)  is  defined  (2)  as 

uo<r)=  Jo(*r;Vo(°fRz)  -  y.<«t)J.(<*rz).  (29) 

Some  of  the  properties  of  this  function  are  discussed  in 
Appendix  C.   Using  Eq.  (29),  Eq.  (27)  becomes 

<hzLr^)=   B*U(*r)SinK«H->)  (30) 

where 

B*  =  B/Ye(<*R2). 

When  boundary  condition  four,  Eq,  (8),  is  applied 
to  Eqs.  (21;)  and  (30),  it  is  seen  that 

AI.(£f?1)S.nKA(H-j.)  =  B*U0(<*R1)5«nhY(H-jJ       (3D 

or  writing  Eq,  (31)  in  a  different  form 

J^SmnACH-^J  -  K2  S^K  *(H-j,).  (32) 

Rearrangement  of  Eq.  (32)  gives 

K  m   5,nKA(H7)  (33) 

Eq.  (33)  implies  that,  regardless  of  the  value  of  y   , 
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the  left-hand  side  of  the  equation  is  always  equal  to 
the  same  constant.  In  order  for  this  to  be  true,  X    and 
X   would  have  to  be  equal.   Thus,  Eq.  (26)  may  be  writ- 
ten as 


(310 


-  «*   +  A2  -  ^  =  O  , 

When  boundary  conditions  four  and  five,  Eqs.  (8) 
and  (9),  are  applied  to  Eqs.  (24)  and  (30) 

AI0(ySR1)  -   B*LUo<Rl)  (35) 

and 

DiA/SI^RD  -  -D2  6*0(4 (*R1).  (36) 

Dividing  Eq.    (36)   by   (35) 

I.C/5R1)  IU<*R1)    .   '  (37) 

Solving  Eq.    (37)    for    /3      gives 

/0  IU«R1)    D,  X,  (ySRD  l3  ' 

Examination  of  Eq.    (38)   reveals  that  it  is  transcen- 
dental in  nature   and  requires   an  iterative  type   of  solu- 
tion to  determine  the  constants    o(     and   /$      .     But  in 
solving  for    o<     and    /3      ,   it  must  be  remembered  that  they 
are  related  by  Eqs.    (16)    and   (3k) i   !•••» 

«* '- X  -  Xl  -(?  .  (39) 
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Because  of  the  periodicity  of  the  functions  J0(r)  and 
Y0(r)  that  compose  the  function  TJ0(r),  an  infinity  of 
solutions  is  obtained  for  Eq.  (38).  Therefore,  Eq.  (39) 
now  becomes 

<  =  <  -  tf  ~/3l  {k0) 

and  the  expressions  for  the  fluxes  are  now  written  as 

OO 

^(r^)=Z!  AnI.(pnr)Si«KAn(H-j,)  (W) 

and 

The  subscript  "n"  on  A  is  because  of  the  infinity  of 
values  of  o<     and  A 

The  application  of  boundary  conditions  six  and 

seven,  Eqs.  (10)  and  (11),  yield 

go 

S  =  D^Av>AY,I0(/5^)CoshAv,M  (U3) 

Ylso 

and 

5  =  D2^B^AnUe^nr)C0sV,  AnH.  (kh) 


n*o 


Utilizing  the  orthogonality  relations  for  Bessel  func- 
tions, the  coefficients  for  Eq.  (k3)    are  obtained  in 


the  usual  manner.  Thus 

oo      /-R1 


T  I  5r  Jo(/SPr;dr  = 

fa]  '  (k$) 

\A^X  A*r  J>(Prr)Ijpnr)C»*\,  A^H  cir 


p»o  neo/_ 
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so  that 

a    -__, SJLCgrLR1)  ,IM 

n"   CjAJJ0(/3.R1)I4(y3wR,1)^J4(^R1)IQ(fi.R1))Cos^AyvH     m) 
z 

if  S  /  S(r).   Similarly,  applying  the  orthogonality 

relations  for  U  (  o(n  r)  as  developed  in  Appendix  C,  the 

coefficients  B  can  be  determined.  Thus 
n 

>   SrU0(*Pr)dr  = 

p-°jRt     -  -   R2  (U7) 

/D*ZZ  /  BtAnrU.C«>r;a(«„r)drJCo,KAnH 

p.o  r»«oJRA 

so  that 

2 

if  3  /  S(r)  and  as  before 

B*  =  Bn/V.(«„R2)  (U9) 

where 

J,  =  -_^rJ4if^Ril  (See  Appendix  c).    (50) 

Another  set  of  source  boundary  conditions  that 
might  be  of  interest  is  one  where  the  source  is  con- 
sidered to  be  zero  at  the  x-y  plane  in  the  absorber, 
while  the  source  in  this  position  in  the  moderator 
remains  as  described  by  Eq.  (11). 

To  describe  the  boundary  conditions  for  this  new 
source  condition,  the  partial  neutron  current  equations 
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are  used.  That  is, 

J+(  <v)    ■  0     for  0  4  r<  Rl   (Region  I)     {$!) 

J+(  y>  )|  m0  ■  S     f or  Rl  <  r  <  R2   (Region  II)    (52) 

where  the  partial  neutron  current  in  the  positive  direc- 
tion is  defined  as  (6) 

Using  Eqs.  (lj.1)  and  (J4.2)  as  definitions  for  <j> 
and  $£  respectively,  the  following  two  equations  are 


obtained: 
00 


4 


+ 


r»=© 
and 


^4»A„I.(y8nrtCo,kAnH  =0 


(5k) 


y   B.U.(«i,r)SlnV>  AnH  + 


ee 


4 
^B*U.tonDA„CoSKAnH  -S. 

Working  first  with  Eq.  (5U)i  it  is  seen  that 


(55) 


XXlo</3wr)(SmU»H  +2D1Av,CosV.AyiH}  =  0,  (56) 

If  fi     is  real,  I  (/5„r)  is  never  equal  to  zero,  and  a 
trivial  solution  is  obtained  if  AQ  »  0.  Therefore,  the 
bracketed  term  must  equal  zero.  After  transposing,  the 
bracketed  term  may  be  written  as 
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To^h  AnU  --ZQ  An  .  (57) 

Since  the  hyperbolic  tangent  is  always  greater  than  or 
equal  to  zero  if  its  argument  is  positive,  the  possible 
solutions  to  Eq.  (57)  are 

When  An-0  ,  solution  of  Eq.  (18)  yields 

Z  =  d  +Ct^,  (59) 

Application  of  Eq.  (7)  leads  to 

Z=  C/1-^/R)  =  C3(M-^)  (60) 

and  the  flux  equation  for  the  absorber  is  written  as 

+<   -£>,,I.flS»rXH-3.)..  (6D 

Utilizing  the  boundary  condition  at  the  interface, 
Eq.  (8) 

^ABX.^R1)(H-^)-J]BtU.C«wkR1)SmK^(Hr^    (62) 


n 


or 


I- 


Since  the  left-hand  side  is  a  constant,  the  above 


implies  that,  regardless  of  the  value  of  ^  ,  the 


(63) 
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right-hand  side  is  Squal  to  a  constant.   For  this  to  be 
true,  ^  must  ec^al  zero  so  that  the  solution  of  Eq. 
(2)  will  yield  th«  same  result  rs  the  solution  of  Eq. 
(7)  when  A^  is  ..qual  to  zero. 

However,  from  physical  consideration  of  the  system, 
it  is  known  that  the  attenuation  of  the  flux  in  the  z 
direction  is  exponential  and  not  linear  as  these  results 
suggest.  Thus,  one  must  conclude  that  the  solution  A^  =  0 
does  not  well  describe  the  physical  situation  and  is 
invalid. 

Now  consider  the  case  when  An  is  imaginary  and 
let  X^-Xfy,   »  then  Eq.  (57)  becomes 

Tan?,H  =  -20,$,  i  (64) 

An  examination  of  a  plot  of  the  tangent  reveals  that 
an  infinity  of  solutions  to  Eq.  (64)  are  possible. 
Because  the  form  of  the  flux  equation  with  A^=^^  is 

*1   =EA*X°(A>r)^s^  <-<H>fl  (65) 

Yl 

negative  imaginary  values  of  i^  are  necessary  to  obtain 
positive  values  of  ^  .  Hence,  the  argument  of  the 
sine  term  must  be  positive. 

Application  of  the  interface  boundary  condition, 
Eq.  (8),  leads  to  the  condition  that 

y  AttXpCflttRi)   _   stnK  yw(H->)  (66) 
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Again,  It  is  evident  that  the  right-hand  side  is  equal 
to  a  constant  which  implies  that  2fh  is  imaginary,  so 
that  once  more  AA  =  ^  • 

The  sine  distribution  in  the  z  direction  does  not 
contradict  the  physical  situation.   Since  the  model  chosen 
calls  for  a  plane  source  of  monoenergetic  neutrons  at  the 
source  plane,  the  value  of  the  neutron  flux  at  the  source 
plane  should  be  largest  and  decrease  to  zero  at  the 
extrapolated  height  of  the  system.   The  sine  function  as 
defined  fulfills  these  conditions. 

To  determine  the  form  of  the  A^  coefficients,  the 
interface  boundary  condition,  Eq.  (8),  must  be  applied, 
and  the  An's  can  then  be  expressed  in  terms  of  the  I^'s. 

Therefore,  the  forms  of  the  coefficients  13  must  be 

*  n 

determined.   Working  with  Eq.  (55)»  rearrangement  yields, 
after  the  substitution  of  sine  and  cosine  for  hyperbolic 
sine  and  cosine, 

4S  =£B*  U.(«nD[Sm  A^H  +  2D2  A^  Cos  XnW]      (67) 

where  the  bracketed  term  is  some  constant,  C  say. 
Applying  the  orthogonality  relation  for  U  (  o<n  r)  as 
developed  in  Appendix  0 

R2 


4sT  rU0(^r)dr« 

£  (RB^  r  U.C*yr)  Uo(*nr0  Cndr 


(68) 


which  beoomes 
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(70) 


B*  =  ^  [(#)'-<*;■  *Jtf/{  R1  U^RI)}2]  ■         (69) 
Thus,  it  is  seen  that  the  coefficients  B"*  are 

similar  in  form  to  those  obtained  from  Eq.  (48).   Now 

Eq.  (8)  is  applied  to  find  AQ: 

£b*U(*wR1)Cc'  Sm  £,<H-j,» 

n 
and 

where  b£  is  defined  by  Eq.  (69)'. 

A  second  approach  to  formulating  this  source  con- 
dition is  by  use  of  the  net  current  equation  at  the 
source  plane;  i.e., 

J(y)\    _o=  O  O   «  r  <  R1 

>_°  (72) 

3<y)\y*o  -  s  R1<   r  <  R2. 

where  J(v)   is  defined   (6)   as 

J(^)  -    -OV  +  ,  r  (73) 

Using  Eq.    (l+l)    and   (2+2)    as  definitions   of    <^    and 
5%      the  following  equations  are  obtained: 


D^Xl.fA,^  A„CosK  AnH=  O 


(74) 
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and 

D27  B*U0(c*„r\)  A*  CosK  AnH  -S  ,  (75) 


Examination  of  Eq.  (74)  reveals  that  trivial  solutions 
are  obtained  if  either  An  =  0  or  A      =  0.   Assuming  that 
S^    is  real,  IQ(  A^   r)  is  never  equal  to  zero.  Thus, 

A  must  be  imaginary  for  a  solution  to  exist.  Letting 

Aw  =  jl  fn  ,   Eq.  (7k)   becomes 

D£A»x°(P»r)uf»Cosf«M   m°"  (76) 

The  roots  of  this  equation  occur  at 

5„  =  (Zk  +  \)-w/ZH  .  (77) 

By  arguments  similar  to  those  used  to  obtain  Eq.  (66), 
it  can  be  shown  that  #y,  must  also  be  imaginary.  Thus, 
Eq.  (75)  becomes 

D2£B*Uo(«n™<;$„Cos  t*H)  =  S,  (78) 

V\ 

Applying  the  orthogonality  condition  for  U0(  of  r) 
as  developed  in  Appendix  C 

■HZ 


s[rlUo<.r)dr  - 

D2V  S'^Cos^H    rU0(*«r)U.(«nr)dr 


(79) 


which  becomes 
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Application  of  the  interface  boundary  condition,  Eq.  (8), 
yields  an  expression  for  An.  Thus, 

*  (81) 

n 

and 

Av^      IeC/3„R1)   "  (82) 

where  B*  is  defined  by  Eq.  (80). 

The  net  current  source  condition  may  be  a  more  accu- 
rate description  of  the  events  that  occur  at  the  source  plane 
than  the  partial  current  source  condition  since  the  net  cur- 
rent equation  accounts  for  neutrons  that  enter  into  the  source 
region  and  then  reenter  the  rod-moderator  region.  The  partial 
current  equation  considers  only  those  neutrons  that  are  moving 

in  the  proper  direction  and  assumes  that  all  events  that 
direct  neutrons  in  the  opposite  direction  are  losses.  The 
fact  that  the  net  current  condition  has  implied  that  an  imag- 
inary source  is  necessary  definitely  will  hinder  its  physical 
application  although  it  is  an  interesting  mathematical  cri- 
terion. Therefore,  the  partial  current  boundary  condition 
would  normally  be  utilized. 

Thus,  for  the  case  of  a  known  plane  isotropic  source  of 
neutrons  located  in  the  x-y  plane,  the  value  of  the  effective 
rod  radius  can  be  determined  analytically  by  solving,  using  a 
trial  and  error  iterative  procedure,  Eqs.  (38),  (ij.0),  (1*6), 
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and  (lj.8).   If  the  value  of  the  source  is  unknown,  the  effec- 
tive rod  radius  can  be  determir  i  numerically  by  fitting  an 
experimentally  determined  flux  plot  by  a  method  of  least 
squares  while  varying  the  rod  radl  l  to  reduce  the  least 
square  error  to  a  minimum.   An  IBM  1620  computer  program  that 
determines  the  effective  rod  radius  and  the  unknown  coeffi- 
cients in  the  flux  equations  is  described  in  Appendix  F. 
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EXPERIMENTAL  PROCEDURE 

To  determine  the  effective  rod  radius,  one  must  first 
obtain,  by  a  suitable  method,  an  experimental  flux  plot  In 
the  absorber  and  moderator  region,  A  method  of  doing  this, 
and  the  method  employed  in  this  work,  Involved  the  use  of 
gold  wire  and  gold  foil  in  conjunction  with  absorber  rods 
of  two  different  diameters  utilized  with  a  graphite  thermal 
column  of  a  nuclear  reactor. 

The  absorbing  rod  used  to  obtain  the  data  for  the  flux 
plot  was  cut  into  three  sections.   Ideally,  these  sections 
were  of  such  a  length  that  the  end  effects  due  to  the  length 
of  the  sections  were  negligible,  and  the  center  section  of 
the  rod  was  of  sufficient  length  that  the  z  dependence  of  the 
flux  was  eliminated. 

Slots  were  cut  into  each  end  of  the  center  section  of 
the  rod,  and  these  slots  were  of  sufficient  depth  and  width 
that  they  accommodated  both  the  foil  and  the  wire.  The  wire 
or  foil  when  in  position  extended  from  the  edge  of  the  spe- 
cial graphite  stringer  completely  through  the  absorber  rod. 
With  the  wire  and  rod  in  place  in  the  stringer,  it  was  then 
placed  into  the  graphite  thermal  column  and  the  irradiation 
period  was  begun. 

For  the  case  where  two  different  rod  diameters  were 
used,  the  stringer  was  constructed  so  that  when  the  one  rod 
was  being  used  to  amass  the  data,  a  graphite  rod  of  the  same 
material  as  the  stringer  was  placed  in  the  groove  for  the 
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second  rod.   The  purpose  of  the  graphite  rod  was  to  eliminate 
the  void  space  that  would  have  been  present  when  only  one 
rod  was  placed  in  the  stringer. 

At  the  conclusion  of  the  period  of  irradiation,  the 
reactor  was  shut  down,  and  the  special  stringer  was  removed 
from  the  thermal  column.  The  stringer  was  then  placed  in  a 
guarded  area  to  permit  the  short  lived  isotopes  of  the  absorber 
rod  to  decay  away.   After  a  suitable  cooling  off  period,  the 
wires  or  foils  were  removed  from  the  rod-stringer  assembly 
and  transported  to  the  counting  facility. 

To  count  the  wire,  a  wire  counting  facility  was  needed. 
This  facility  consisted  of  a  wire  guide  to  measure  the  posi- 
tion of  the  wire  in  relation  to  the  center  of  the  rod  and 
a  Geiger  counting  tube.  The  counting  tube  was  shielded  to 
prevent  the  beta  radiation  from  scattering  and  streaming  into 
the  sensitive  volume  of  the  tube. 

A  foil  counting  facility  similar  to  the  one  constructed 
to  count  the  wire  was  feasible,  but  a  more  simple  arrangement 
was  employed.  The  foil  was  taped  to  a  piece  of  graph  paper 
and  then  cut  into  small  strips  at  appropriate  intervals. 
These  strips  were  then  placed  on  a  special  plancet  that  insured 
reproducible  geometry  with  respect  to  the  counting  tube,  and 
the  small  foils  were  then  counted  using  a  gas  flow  proportional 
counter. 
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ANALYSIS  OP  DATA 

In  the  theoretical  development  of  the  two  region  P-, 
transport  model,  it  was  shown  that  an  analytic  expression 
for  the  neutron  flux  in  both  regions  can  be  determined  for  a 
given  absorber  radius  once  a  source  condition  is  established. 
However,  to  determine  the  effective  rod  radius,  it  was  nec- 
essary to  fit  a  set  of  experimental  data  with  the  proper 
expressions  for  the  flux  while  varying  the  absorber  radius. 
An  IBM  1620  computer  program  was  written  for  this  purpose, 
and  it  is  described  in  Appendix  P. 

The  effective  rod  radius  was  determined  by  a  trial  and 
error  procedure  which  strived  to  minimize  the  sura  of  the 
weighted  squares  of  the  residuals  between  the  experimental 
values  and  the  values  calculated  from  the  theoretical  expres- 
sions for  the  neutron  flux  in  both  regions.  The  weighting 
function  used  was  that  one  which  required  that  the  residual 
between  the  calculated  value  at  any  point  and  the  experimental 
measurement  should  lie  within  the  statistical  standard  devia- 
tion. 

Another  point  considered  in  the  fitting  of  the  experi- 
mental data  was  the  number  of  harmonics  necessary  to  best 
fit  the  data.   Provision  was  made  in  the  computer  program  to 
accommodate  up  to  five  harmonics.  However,  the  inclusion  of 
more  than  three  harmonics  was  deemed  inadvisable  since  the 
accuracy  of  the  data  did  not  warrant  a  more  precise  analysis. 

All  of  the  data  were  analyzed  using  one  harmonic,  and 
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two  different  situations  were  considered  in  the  analysis.   In 
one  case,  the  value  obtained  from  the  solution  of  the  trans- 
cendental equation  <K/X  =  tanh  M/Zs  was  used  for  X     in  the 

rod;  in  the  second  case,  the  diffusion  theory  value  for  X  in 

i 
the  rod,  X     =  (3>1l<^T.5    )    $   was  employed.   The  purpose  of  this 

approach  was  to  show  the  difference  between  asymptotic  theory 

and  the  P-j_  transport  theory  in  the  absorbing  media,  since 

diffusion  theory  is  true  only  when  T.^  «   Es 

Three  particular  sets  of  data  were  analyzed  using  two 

and  three  harmonics  to  see  what  effect,  if  any,  the  additional 

terms  would  have  on  the  value  of  the  least  square  error.   In 

each  of  these  cases,  the  value  used  for  }£  in  the  rod  was  that 

obtained  from  asymptotic  theory. 
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RESULTS  AND  CONCLUSIONS 

The  analysis  of  the  data  for  this  work  was  identified 
by  the  following  notation  and  parallels  the  notation  used  in 
the  tables  of  original  data  found  in  Appendix  D: 

Case  1.   One  inch  rod  with  gold  wire  near  source. 

Case  2.   One  inch  rod  with  gold  wire  away  from  source. 

Case  3.   One  inch  rod  with  gold  foil  near  source. 

Case  4.   One  inch  rod  with  gold  foil  away  from  source. 

Case  5.   Two  inch  rod  with  gold  wire  near  source. 

Case  6.   Two  inch  rod  with  gold  wire  away  from  source. 
The  phrases  "near  source"  and  "away  from  source"  refer  to 
the  position  of  the  gold  foil  or  wire  with  respect  to  the 
reactor  core  during  the  period  of  irradiation.   Thus,  "near 
source"  signifies  that  the  foil  (or  wire)  was  placed  in  the 
special  graphite  stringer  in  the  position  nearest  the  core, 
and  "away  from  source"  signifies  that  the  foil  (or  wire)  was 
placed  in  the  special  graphite  stringer  in  the  position  far- 
thest from  the  source  (See  Figure  23  in  Appendix  D). 

The  Tabulation  of  Results  was  further  subdivided  by  the 
letters  A,  B,  C,  and  D  to  indicate  the  method  of  analysis. 
These  letters  have  the  following  meaning: 

A.  Data  analyzed  using  one  harmonic  and  the  value  of 
Ki  obtained  from  asymptotic  theory. 

B.  Data  analyzed  using  one  harmonic  and  the  value  of 
Ki  obtained  from  diffusion  theory. 

C.  Data  analyzed  using  two  harmonics  and  the  value  of 
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^   obtained  from  asymptotic  theory, 

D.   Data  analyzed  using  three  harmonics  and  the  value 
of  XA  obtained  from  asymptotic  theory. 

For  example,  the  results  listed  as  Case  3B  signify  that 
the  data  recorded  using  the  gold  foil  positioned  in  the  stringer 
nearest  the  core  were  analyzed  using  one  harmonic  and  the 
diffusion  theory  value  for  X   in  the  rod. 

Table  1  presents  a  tabulation  of  the  results,  identified 
by  the  preceding  notation,  obtained  for  the  value  of  the 
effective  rod  radius,  Reff  >  for  both  the  one  inch  and  the  two 
inch  diameter  steel  rod  using  the  IBM  1620  computer  program 
described  in  Appendix  F.   Figures  2  through  20  graphically 
illustrate  the  variation  of  the  flux  in  the  absorbing  rod 
and  moderator  calculated  by  the  least  square  analysis  as  com- 
pared to  the  experimental  data. 

For  the  analysis  that  utilized  the  value  of  X   in  the 
rod  obtained  from  the  solution  of  the  transcendental  equation 
X/Z  =  tanh  X/Zs,  in  only  two  instances  did  the  value  obtained 
for  Rgff  exceed  the  value  of  the  physical  rod  radius.   Both  of 
these  cases  involved  the  one  inch  rod.   In  one  Instance,  spe- 
cifically Case  4,  the  value  obtained  for  R  ff   was  not  consist- 
ent with  the  values  obtained  in  the  other  cases  for  the  one 
inch  rod.   However,  the  validity  of  this  data  was  questioned 
at  the  time  it  was  recorded  because  the  gold  foil  was  noted 
to  be  out  of  position  when  it  was  removed  from  the  special 
stringer.   In  both  cases  involving  the  two  Inch  rod,  the  value 
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of  Ref;C  was  less  than  the  rod  rcdius.  Thus,  one  must  conclude 
from  these  results  that  this  method  of  analysis  was  not  too 
satisfying  from  the  stand  point  of  an  effective  rod  radius. 
Although  the  calculated  rod  radius  was  greater  than  the  phys- 
ical radius  of  the  rod  in  two  instances,  the  majority  of  the 
results  led  to  the  conclusion  that  the  absorbing  rod  was  smaller 
in  size  than  its  actual  measurement.  Thus,  the  accuracy  of 
the  experimental  data  has  left  something  to  be  desired. 

In  the  analysis  that  used  the  diffusion  theory  value  for 
X  in  the  rod,  the  value  obtained  for  Reff  in  all  cases  was 
smaller  than  the  rod  radius.  This  result  was  not  too  surpris- 
ing since  diffusion  theory  is  known  (11)  to  be  faulty  in 
regions  where  ZT<^  is  of  the  same  order  of  magnitude  as  I  . 
Hence,  this  method  of  calculating  the  effective  rod  radius 
was  less  satisfying  than  the  first  and  should  not  be  used  to 
determine  the  effective  rod  radius. 

When  more  than  one  harmonic  was  included  in  the  analysis 
of  the  data,  the  general  result  was  that  the  value  of  Reff 
increased  and  the  least  square  error  decreased.  In  one  par- 
ticular case  utilizing  the  one  inch  rod,  Case  1,  the  value 
of  Reff  calculated  using  one  harmonic  increased  from  0. lj.60lp 
inch  to  0.5063  inch  by  using  three  harmonics  while  the  least 
square  error  decreased  from  5^-1  *k   "to  352.1  respectively. 
The  results  of  Case  2  present  a  more  definite  illustration 
of  how  the  additional  terms  of  the  multi-harmonic  analysis 
reduced  the  value  of  the  least  square  error.  For  this  case, 

the  value  of  R  ~~   calculated  using  one  harmonic  increased 
eff 
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from  0.5072  inch  to  0.5^88  inch  for  two  harmonics  and  then 
decreased  to  O.5138  inch  for  three  harmonics.  The  value  of 
the  least  square  error  in  these  three  instances  varied  from 
23i|.8  to  88.5  to  99.8  respectively. 

The  result  that  the  inclusion  of  two  or  three  harmonics 
improved  the  value  of  the  least  square  error  indicates  that 
the  foil  or  wire  were  too  near  the  source  and/or  the  end  of 
the  rod,  and  this  resulted  in  the  significance  of  the  higher 
order  harmonic  terms.   If  the  foil  or  wire  had  been  located 
at  a  greater  distance  from  the  reactor  core  and  if  the  rod 
had  been  longer,  these  terms  would  have  been  insignificant 
compared  to  the  fundamental  mode,  one  harmonic. 

Case  2  results  emphasize  this  fact.  The  wire  used  to 
record  the  data  was  located  away  from  the  core  which  means 
that  it  was  ten  inches  farther  from  the  core  than  was  the 
wire  that  was  located  near  the  core,  Case  1.  For  Case  2  the 
minimum  least  square  error  and  the  maximum  value  for  RQff 
occurred  using  two  harmonics  to  analyze  the  data.  Thus,  it 
appears  that  this  little  increase  in  distance  was  sufficient 
to  reduce  the  effect  of  the  third  harmonic  on  the  value  of 
the  flux. 

The  results  for  the  multi-harmonic  analysis  of  the  two 
inch  rod  data  were  totally  unexpected.   Case  6  was  selected 
for  analysis  since  it  presented  the  best  prospect  for  pro- 
ducing an  RQff.  on  the  order  of  the  physical  dimension  of  the 
rod  judging  from  the  results  obtained  for  Case  1  and  Case  2. 
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For  two  harmonics,  the  desired  result  was  obtained  when  the 

value  of  the  least  square  error  decreased  and  R  ~~   increased 

eri 

from  the  values  obtained  for  one  harmonic.  However,  when  the 
three  harmonic  analysis  was  attempted,  the  value  of  both  the 
least  square  error  and  RQff  decreased.  Thus,  it  appears  that, 
as  with  Case  2,  the  wire  was  located  far  enough  from  the  core 
to  reduce  the  harmonic  content  of  the  flux  to  two  harmonics, 
but  the  data  were  insufficiently  accurate  to  produce  the 
expected  results. 

It  was  also  noted  in  the  results  of  the  multi-harmonic 
analysis  that  the  sign  on  the  term  associated  with  the  fun- 
damental mode  was  negative  in  all  cases,  and  that  the  second 
harmonic  was  the  dominant  term.  Referring  to  Appendix  G,  it 
can  be  seen  that  the  signs  on  the  coefficients  of  the  flux 
equations  are  exactly  opposite  what  one's  intuition  would 
suggest;  i.e.,  one  would  expect  the  fundamental  mode,  the 
first  harmonic,  to  be  dominant  and  have  positive  coefficients, 
while  the  higher  harmonics  would  be  smaller  in  magnitude  and 
have  either  positive  or  negative  signs.  Hence,  it  would  appear 
that  the  original  data  were  not  as  accurate  as  they  could  have 
been. 

In  the  graphs  of  the  original  data,  it  was  noticed  from 
the  beginning  that  unexpected  peaks  were  evident  in  the  moder- 
ating region.  At  first  it  was  thought  that  these  peaks  were 
due  to  the  location  of  the  control  rods  used  to  bring  the 
reactor  to  critical.  But,  in  a  later  irradiation,  these  rods 


were  withdrawn  completely  and  different  ones  were  employed  to 
cause  the  reactor  to  become  critical  and  the  peaks  were  still 
present.   No  explanation  for  these  peaks  was  found  when  the  data 
were  recorded,  and  none  is  off  ere.,  at  this  time. 

Summarizing  the  results,  it  was  seen  that  values  for 
Reff  greater  than  the  physical  dimension  of  the  rod  were  ob- 
tained in  a  few  cases,  and  in  these  cases,  the  magnitude  of 
the  flux  depression  was  very  close  to  that  obtained  experi- 
mentally. For  the  cases  in  which  Ref^  was  found  to  be  less 
than  the  actual  rod  dimension,  in  one  instance  the  data  were 
known  to  be  invalid,  and  in  the  others,  the  results  were 
inconclusive  since  a  multi-harmonic  analysis  of  the  data  was 

not  conducted. 

Thus,  it  would  appear  that  the  concept  of  an  effective 
rod  radius  was  reasonable,  but  it  will  require  much  more 
investigation  to  produce  graphs  of  the  type  that  are  needed 
to  fully  utilize  this  concept.  Great  care  must  be  exercised 
when  data  for  this  type  of  experiment  are  recorded,  and  the 
location  of  the  detecting  material  in  relation  to  the  reactor 
core  is  also  important.  To  analyze  the  data  for  this  type  of 
investigation,  computers  of  greater  speed  than  the  IBM  1620 
would  be  very  advantageous,  especially  for  the  multi-harmonic 
analysis.   Examination  of  other  materials  and  various  rod 
sizes  using  the  ideas  embodied  in  this  work  would  provide  an 
excellent  starting  place  for  future  experimentation. 
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Fig. 3      Flux    Calculated  by  Least  Squares  Analysis 
Case    IB 
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Fig. 4      Flux    Calculated   by  Least  Squares   Analysis 
Case    IC 
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Fig.  5      Flux   Calculated   by  Least  Squares  Analysis 
Case   ID 
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Fig. 7       Flux    Calculated    by   Least   Squares   Analysis 
Case    2B 
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Fig.  8        Flux  Calculated  by   Least    Squares    Analysis 
Case    2C 
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Fig. 9       Flux   Calculated  by    Least    Squares    Analysis 
Case     2D 
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Fig.ll       Flux    Calculated    by    Least   Squares  Analysis 
Case     3B 
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G     Experimental    Points 
Effective    Rod  Radius 


0  1  J. 

4  2 

Distance    from  center    of   one    inch  rod   (inches) 

Fig.20    Comparison    of    Flux   in   Rod    Calculated  by 
Least    Squares    Analysis  for    Case    I 
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Suggestions  for  Further  Work 

For  the  concept  of  an  effective  rod  radius  to  be  of 
use,  graphs  relating  the  dimensionless  parameter  Z!^R1  to 
the  effective  rod  radius  should  be  available.  To  construct 
graphs  of  this  type  and  to  determine  the  empirical  represen- 
tation for  these  graphs,  various  combinations  of  rod  radius 
and  absorber  material  should  be  investigated. 

Since  the  primary  problem  in  the  determination  of  the 
effective  rod  radius  is  in  the  calculation  of  the  coeffi- 
cients in  the  flux  equations,  other  representations  for  these 
terms  could  be  examined.  More  rigorous  boundary  conditions, 
e.g.,  no  return  current  at  the  outer  boundary,  might  be  exam- 
ined to  see  what  effect,  if  any,  they  would  have  on  the  value 
of  the  effective  rod  radius. 

Because  P^  transport  theory  is  just  slightly  more  accu- 
rate than  elementary  diffusion  theory,  investigations  consid- 
ering higher  order  P  calculations  should  prove  of  interest. 
Some  work  has  already  been  done  in  this  area  (13) ,  but  there 
is  still  much  to  be  accomplished.  This  type  of  analysis  is 
somewhat  more  difficult,  but  the  improvement  in  the  accuracy 
of  the  results  should  more  than  compensate  for  the  extra 
labor  involved. 
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APPENDIX  A 

Derivation  of  the  Boltzmann  Transport  Equation 
in  Cylindrical  Geometry 

Following  the  notation  of  Weinberg  and  Wigner  (16), 
the  Boltzmann  transport  equation  for  the  monoenergetic, 
steady  state  system  is  given  as 


(A-l) 


where 

fl    is  the  unit  vector  in  the  direction  of  the 

neutron  motion, 
f(t,D,)  drd A    is  the  angular  flux;  i.e.,  it  is 
the  number  of  neutrons  in  the  volume  element 
dr  around  r  whose  directions  of  motion  lie  in 
the  solid  angle  6Q.   about  X2  times  the  speed  v, 
dH   is  the  differential  solid  angle, 
Z    is  the  total  cross  section  (Z  =  £<*.  +  Zs  )     , 
£s(il~'>'Q)     is  the  scattering  cross  section  for 

changing  XV   into  the  direction  dfl   about  £L  , 
Q-Qrad  -f  is  the  directional  derivative  of  f  in 
the  direction  £X    ;   i.e.,  if  the  coordinate  s 
is  extended  in  the  direction  =0.   ,  then 

r\      ~  j  £  =  ii  <A-2> 
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For  a  system  with  cylindrical  symmetry  and  which  is 
infinite  in  the  z  direction,  the  angular  flux,  f,  depends 
on  three  variables.   In  the  notation  of  Weinberg  and 
Wigner  (16),  these  variables  are: 

r  —  the  distance  from  the  axis  of  the  cylinder, 

the  z  axis 
Q   --  the  angle  between  £l    and  z;  and 
<f   --  the  angle  between  the  projection  of  =Q=  on 
the  x-y  plane  and  the  direction  of  r. 
The  geometry  of  the  above  described  system  is  shown  in 
Figure  21. 


ds  Cos  0 


?+<}<f> 


Figure  21.   Cylindrical  Geometry 
Coordinate  System 


Note  that  when  d  ¥     is  small 


5m(P+<iW  =    Sm  9 


Qos(9>*6P)   =    Cos  V  , 
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Thus, 

-r  ^  =  as  5m  e  Sm  7  (a-3) 

and 

dr  =   ds  Sm  9  CosV  ,  (A-U) 

Then 

7^  =   -~   Smd  SiwV  (A-5) 

and 

^£  =    Sm©  Cos9>,  (A-6) 


v  r 
Prom  Eq.    (A-2)        Q  •  qrad  -P  =  -r—    - 

Since 

f  =  f(r,  9>)  (a-7) 

then 

ds       ^r    ds        ay    ds  (A     J 

so  that 

HBi£[Sin©Cb.^]  ^4ll--rs^a5'^]  (A~9) 

^ s      <^rL  ^  L   r  ' 

The  substitution  of  Eq.  (A-8)  and  Eq.  (A-9)  into 
Eq.  (A-l)  yields 


S 
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(A-10) 


/< 


djQ  Zs(r,cos9jfCr;^^ 

which  is  the  monoenergetic  Boltzmann  transport  equation 
in  cylindrical  geometry  with  y     and  ¥     symmetry. 

To  transform  the  Boltzmann  equation  to  a  more  con- 
venient form  for  computation,  expand  Zs(r;  Cos^o),  $(r9,¥) 
and  S(T}  0,  t)    in  spherical  harmonics  noting  that  Hs 
only  depends  on  Qa     ,  the  angle  between  XI  and  Q,' 
Thus 

«i;e,w-£^  ^?J^>  <A"12) 

s^e.wtt  &'?w.<-to  (A"13) 

where  one  form  of  P  U20  is  given  by 

In  order  that  the  spherical  harmonics  can  be  fully 
utilized,  the  following  relations  for  this  special  use 
should  be  noted: 

P.     (Ik)   =<-Dlr"i?M  (H)*  (A-1*) 
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where  the   indicates  the  complex  conjugate; 

the  orthogonality  relation  for  Legendre  polynomials; 

jdfl^ifD^iSLf  -2^£/.*»w  <A"l6> 

the  addition  theorem  for  Legendre  polynomials; 

X 

5<Cca>-£g„ui>5rf<ijQ*    ■  (A-i7) 

and  finally,  the  recursion  relations  for  the  Legendre 
polynomials 

.    s •"  ee<v -  ^^(^-1)9.?-^  -2ffl-'*aM6liWft)l      (A"18) 

where 


i 


^"  £^*3)  ^-  2f2i  +  3) 


The  substitution  of  Eqs.  (A-ll),  (A-12),  and  (A-13) 
into  Eq.  (A-10)  leads  to 

OO      A  / 


EI  /i^SLUDZfr)  =££  S^tli)  (A-20) 


Now  multiply  Eq.  (A-20)  by  R.Ui)  dfl   aod  integrate 
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over  all  Q   .   Examining  each  term,  starting  with  the 
right-hand  side,  keeping  in  mind  the  orthogonality 
relation,  it  can  be  shown  that: 


[47^(ti+l8f  tr)< 


>(. 


JLL  s^uui^a^-Wrw+uls^  (1-22) 


,.  00    X 


LL  W$„WZM%J&>  i&  •forfu4fJlwto .       (i-23) 

For  the  first  term  on  Eq.  (A-20),  write  Cos?   in  expo- 
nential form  and  then  substitute  the  recursion  relations 
Eqs.  (A-18)  and  (A-19).   Thus 


The  second  term  of  Eq..  (A-20)  is  handled  in  a  like  man- 
ner. 
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Prom  Eq.  (A-14)  it  is  seen  that 

Jp~  -  ^Owg^Ul).  (A-25) 

Thus 


1  ZSJ& 


J 


/J-  smOBe^^ytaZI  t,wa..«t)*dil-i*.  5™(^  = 


Regrouping  Eqs.  (A-21),  (A-22),  (A-24)  and  (A-25)  in  the 
proper  order,  the  spherical  harmonics  form  for  the 
Boltzmann  transport  equation  in  cylindrical  coordi- 
nates is  obtained 
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APPENDIX  B 


Derivation  of  the  P-^  Transport  Equation 
in  Cylindrical  Geometry 


In  order  to  derive  the  P^  transport  model  for 
cylindrical  geometry,  the  Boltzmann  transport  equation 
must  be  utilized.  The  solution  of  this  equation  can 
be  found  in  many  references  (3,  5*  7»  11*  12,  16).  The 
form  of  the  Boltzmann  equation  that  will  be  used  in  the 
derivation  is  that  one  found  in  Weinberg  and  Wigner  (16) 
page  276  for  a  cylinder  infinite  in  the  z  direction  and 
for  monoenergetic  neutrons.  The  equation  has  the  form 

a     r^    .  mrvH-i  p  _  R    r^ vn-iip  _p    r  ^     ,  gitJl  /> 

(B-l) 


where 


_  /uw)(jwd},/2  R      {U-^+DU-Yn+Zii11 


A^- 272773) <*»        "2TOT37 

**»' — ITzTT) Im~       mi-i) 

For  the  P-j_  approximation,  J.-09^      while  -J.  <7n<+JL 
The  £    equals  zero  for  J.>,Z    ,  and  when  jt  =  —  4 
Because  the  medium  of  interest  is  assumed  to  be  purely 
diffusing,  the  source  term,  S^w   ,  equals  zero.  In 
the  Pn  approximation,  scattering  is  assumed  to  be  lin- 
early anisotropic  so  that  sQ  f   0  and  s^  ^   0  while  all 
other  Sg  =  0 
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Because  of  the  cylindrical  symmetry  around  the 
z  axis,  the  only  £   of  interest  are  those  where 
Jl   +777=  even.   All  other  -A.   are  neglected.   (The 
equations  for  Jl  +  M  -   even  are  independent  of  the 
^  +  Yn  =   odd  so  that  no  part  of  a  particular  solution 
is  omitted. ) 

Keeping  the  preceding  facts  in  mind  along  with  the 
property  that,  because  of  symmetry  in  the  $*     direction, 
fi  *-  a  ("1)  "fi    »  the  following  equations  are  obtained: 

JL-m  Jcyn 

When  2=0,      y„  =  0 

?f|r  -  fKu   +   4^s.-f..  -  Z*.   =  0   .  (B-2) 

When  2.   =  It  >n  =  1 

f#^^^-Zfn=0.  (B-3) 

Solving  Eq.    (B-3)  for  fxl 

3 

and  using  the  definition  for  Z.   that  is  found  in  (16), 

Eq.  (B-4)  becomes 

W    &foo    1 

•t 


£  -  1L!  Use  _L 


Substitution  of  Eq.  (B-5)  into  Eq.  (B-2)  and  assuming 
Z!+r  is  spacially  independent  leads  to 

3Z"  ^  +  "F  !rJf«  +  <4-ts. -Z)f..  =  0.     (B-6) 
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Noting  that  4-TT  s0  -  Z  =  -  Z^  and  defining 

Ds3i;  (B-7) 

tr 

for  the   P1  approximation,    Eq.    (B-6)   becomes 

DV2f„  -  E„f„    =  0  (b-8) 

where 

ttA,    -    4-TTf..  (b-9) 

and  V2  is  the  well  known  Laplaoian  operator. 

Thus,  the  P,  transport  equation  is  identical  in 
form  with  the  thermal  diffusion  equation  except  for  the 
definition  of  D,  the  diffusion  coefficient. 

The  net  current  in  the  r  direction  can  also  be 
obtained  using  the  preceding  method  and  equations. 
When  JL  -   1,  yn   =  -1,  Eq.  (B-l)  becomes 

"  *F  *k   +  Hp  *v.  "  "t-i  "  °  ■  (B"10) 

Now  subtract  Eq.  (B-10)  from  Eq.  (B-3)  and  obtain 

Using  the  previous  definition  of  ZL+  ,  Eq.  (B-ll)  becomes 

/r^-=  -  zi^Lf,,-^]  =o.  (B-i2) 

Now  use  Eq.  (B-7)  to  define  D  and  rearrange  Eq.  (B-12) 
into  the  form 

J(r)  =  -  DV€0  -  4~2^m  "fn^  (B-1^) 
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where  V  represents  the  gradient  operator  in  cylin- 
drical coordinates  with  no  cp   or  z  dependence. 
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APPENDIX  C 
The  Orthogonality  Condition  for  UQ(  o<   r) 
U0(  o(   r)  is  a  solution  of  Bessel»s  equation 

f  JF(rJF>  +  «l~  =  °   .  {c-1] 

where  the  boundary  condition 

GL>(b)  =  0  (C-2) 

has  been  applied.  For  this  case 

co(r)=VJc<r)=   J.C<*r;Y«G*b) -Yo(*rJdft(<*b;  .     C&-3) 

Let  U0(o<   r)   be   a   solution   of  Eq.    (C-l)    and  let 
V0(B  r)   be   a  solution  of 

where    ^/-(r,>  =  U0C/3rj  .     Multiply  Eq.    (C-l)   by    \r  and  Eq. 
(C-4)   by   cj    and  subtract  Eq.    (C-l}.)    from  Eq.    (C-l).     Thus 

(c-5) 

Now  multiply  Eq.  (C-5)  hy  rdr  and  integrate  over  (a,b), 
the  interval  of  interest. 
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where 


(C-7) 


Thus 

(C-6) 

From  the  preceding  definitions  of  I/-  and  dc>  it  is 

easily  seen  that  both  I/*  and  £c>  are  zero  when  evaluated 
at  the  upper  limit.  Thus,  evaluating  Eq.  (C-6)  at  the 
lower  limit 

rb 
(/32-^;Jr^codr  = 

-'a. 

»afUo(pa)[-<vUo(0fa)]-Uo(^a)[T8Uo(j3a^} 

where 

and 

jf  =  -/5[J1(/3r)Y0(/5b)-  ^(^r)  J0(/8b)] 

=  -j8Ue(/8r)  =-/?U1(/sr), 

Eq.  (C-7)  can  be  rearranged  into  the  form 

(/3W)  r^oodr  =  AW^Kl  +  ADoW]  (c-io) 


(C-8) 


(C-9) 


Jk0  =   aUo(/Sa)  (C-ll) 

^>n  =  -/3U0'(/3a;/U0(/3a)  (0-12) 
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Now,  if  (X  and  fi      are  positive  roots  of  the  equation 

<*U'(#<0  ■*  JiiUoioUk)   =  0  (0-13) 


then 

rvcodr   =/ rUof/3r)U0Co<r;dr  =0  (o-i4) 

which  is  the  orthogonality  condition  for  U  (  ©<  r). 
In  order  to  evaluate  the  integral 

-'a 

take  Eq.  (0-1)  and  multiply  it  by  £  r  7^  : 

2rJF-JF(r5T)+  «»-jfrqf*«  =0,       (c-15) 

Rearranging  the  terms  of  Eq.  (C-15)  yields 

d?(r  JF)Z    +  *^  jp(^0")  -  0.  (0-16) 

Now  multiply  Eq.    (C-16)   by  dr  and  integrate   over   (a,b), 
the  region  of  orthogonality: 

Ijr(r37)ldr   +  ^/r^Cw'Jdr    =   0.  (0-17) 

Integrating  the  second  term  by  parts  and  the  straight 
forward  integration  of  the  first  term  yields 

b      „_     _.b 


[Cr^)2^b  +  oc*£rW]b  -icfJruj-Ar  =  0  .  (0-18) 


Rearranging  Eq.  (C-18)  and  evaluating  the  expressions 
for  CO7-   and  lj"p;   at  the  limits  yields  the  following 
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equation: 


2cff  r  U0Vr)dr  = 


(C-19) 


From  Murray  (12),  the  Wronskian  relation  is  given  as 

Jo(X)  \{X)  -  J,(X)  Yo(X)  =  -  tj%  .  (c-20) 

Using  Eq.  (C-8)  as  a  definition  for  U^  oCb),  It  is 
easily  seen  that 

^W-^E  (c"21) 

and  by  applying  the  definition  for  U-,  (  &<   a)  as  given 
in  Eq.  (C-13),  Eq.  (C-19)  becomes 

*  (C-22) 

Another  integral  that  must  be   considered   is 

|rU0(o<r)ar  .  (C-23) 

-"a 

Writing  the  equivalent   expression  for  TJ  (  ex  r) ,    the 
integral  becomes 

Ix-U^oirJdr  =/ r[s)#(ecr)YoCofb)-Y0(ofr)  J0(ofb)]dr  .  (c-210 

Performing  the  indicated  integration: 
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frUotaDdr  =(£[j>r)Y0(*l,J  -\(«r)  JaC*b)]]b 

b 


-[^Ui^;l 


(0-25) 


=  ^lV*b)  -  —  u^Ccxa)  , 


Substituting   the  value  for  U   (  o{  b)   from  Eq.    (C-21)  and 
the  value  for  U1(  <X  a)   from  Eq..    (C-13),    Eq.    (C-25) 
becomes 

J  rU0(cxr)dr  =   ^,  +  ^^UoUa)   ,  (0-26) 

Further  discussion  of  the  properties  of  this  function 
can  be  found  in  Oarslaw  and  Jaeger  (2). 


74 


APPENDIX  D 
Experimental  Facilities 

The  experimental  facilities  used  to  compile  the  prelim- 
inary data  employed  in  this  work  were  the  Argonaut  Training 
Reactor  of  the  International  School  of  Nuclear  Science  and 
Engineering  located  at  Argonne  National  Laboratory.   The  use 
of  these  facilities  was  made  possible  through  the  courtesy 
and  cooperation  of  the  International  School  of  Nuclear  Science 
and  Engineering  and  the  Association  of  Midwestern  Universities. 
It  is  noted  that  the  experimental  phase  of  this  work  was 
completed  before  the  Kansas  State  University  TRIGA  MARK  II 
Reactor  went  critical. 

A  complete  description  of  the  Argonaut  reactor  can  be 
found  in  references  (10,  15).   The  following  will  be  a  short 
account  of  the  actual  facilities  used. 

The  foils  were  irradiated  in  the  external  graphite  ther- 
mal column.   The  thermal  column  is  an  extension  of  the  reactor 
reflector  graphite,  and  it  is  four  feet  by  five  feet  in  cross 
section  by  six  feet  long.   It  is  composed  of  assorted  lengths 
of  machined  reactor  grade  graphite  four  inches  in  cross  sec- 
tion.  A  total  of  fifteen  removable  horizontal  stringers  are 
available  for  access  in  the  direction  of  the  core.   The  J-10 
stringer  was  the  one  in  which  the  experimental  measurements 
were  made.   It  is  located  on  the  vertical  and  horizontal 
center  line  of  the  core. 
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The  core  of  the  Argonaut  reactor  system  is  annular  In 
design,  and  it  is  reflected  internally  and  externally  with 
graphite.   The  fuel  elements  are  hetrogeneous  and  consist  of 
BORAX- type  fuel  plates  (l£)  in  rectangular  assemblies  three 
"by  six  inches.   Light  water  serves  as  the  moderator  in  and 
surrounding  each  fuel  box  while  graphite  moderator  pieces 
separate  the  individual  fuel  boxes. 

Because  of  the  annular  design  of  the  core,  the  core  can 
be  loaded  with  fuel  in  two  different  configurations  for  power 
runs:   a  one  slab  loading  or  a  two  slab  loading.  The  one 
slab  loading  consists  of  six  fuel  boxes  placed  in  the  annulus 
separated  only  by  the  graphite  moderator  pieces.  When  a  two 
slab  loading  is  required,  an  identical  configuration  of  fuel 
elements  is  located  diametrically  opposite  the  initial  load- 
ing. 

The  critical  loading  of  the  two  slab  system  is  approxi- 
mately four  kilograms  of  U-235  at  20$  enrichment.  With  this 
loading  the  reactor  will  normally  operate  in  the  range  of  one 
to  one  hundred  watts.  The  reactor  can  operate  continuously 
at  one  kilowatt  and  is  capable  of  operation  at  a  maximum  ther- 
mal power  of  ten  kilowatts. 

For  this  particular  work,  a  one  slab  loading  was  used, 
and  the  fuel  elements  were  placed  in  position  nearest  the 
thermal  column;  as  shown  in  Figure  22. 
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Experimental  Procedure 

The  data  compiled  for  this  work  were  taken  by  the 
method  outlined  in  the  main  body,  and  the  pertinent  details 
will  be  presented  here. 

The  absorber  rods  were  steel  measuring  one  inch  and  two 
inches  in  diameter.   The  gold  wire  was  used  to  record  data 
for  both  the  one  inch  and  the  two  inch  rod  while  the  gold 
foil  was  used  only  in  conjunction  with  the  one  inch  rod. 
Both  rods  were  cut  into  three  sections;  the  center  section 
was  ten  inches  long;  the  forward  section,  the  section  nearest 
the  core,  was  six  inches  long;  and  the  aft  section  was  eight 
inches  long  (see  figure  23). 

The  center  section  had  a  slot,  1/32  inch  wide  by  9/32 
inch  deep,  located  on  the  center  line  of  the  rod  cut  into 
each  end.   Into  each  slot,  one  gold  wire  or  one  gold  foil 
was  placed  depending  on  the  phase  of  the  experiment.   The 
wire  or  foil  was  of  such  a  length  that  it  extended  from  the 
edge  of  the  special  graphite  stringer  to  approximately  one 
inch  beyond  the  edge  of  the  rod  when  the  rod  was  in  place  in 
the  stringer.   The  special  stringer,  with  a  graphite  rod  in 
the  position  of  the  steel  rod  not  being  used,  was  then  placed 
in  position  for  the  irradiation*   The  wire  or  foil  was  irra- 
diated for  a  period  of  one  hour  while  the  reactor  was  oper- 
ated at  a  power  level  of  100  watts. 

The  reactor  was  shut  down  at  the  conclusion  of  the 
irradiation  period  by  a  "scram"  procedure,  and  the  time  at 
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which  this  occurred  was  recorded.  Following  the  shut  down, 
the  special  stringer  was  removed  from  the  thermal  column, 
and  allowed  to  cool  off.  The  wires,  or  foils,  were  then 
removed  from  the  stringer  and  transported  to  the  counting 
facility  which  was  located  at  some  distance  from  the  irradi- 
ation facility. 

The  wire  was  counted  as  described  previously.  The  count- 
ing tube  was  an  end  window  Geiger-Mueller  type.  As  the  wire 
was  advanced  under  the  G-M  tube,  it  was  always  moved  a  suffi- 
cient distance  to  prevent  counting  a  portion  of  the  same  sec- 
tion of  the  wire  twice.  At  large  distances  from  the  center 
of  the  rod,  the  wire  was  advanced  in  steps  of  one-eighth  inch, 
In  the  immediate  vicinity  of  the  rod,  the  step  was  reduced 
to  one-sixteenth  inch,  the  diameter  of  the  hole  in  the  wire 
guide.  This  increment  was  maintained  until  the  wire  was 
advanced  past  the  region  where  the  rod  had  been  positioned. 
At  each  position  where  a  count  was  taken,  the  time  when  the 
counting  began  was  recorded. 

The  foil  counting  procedure  was  that  one  described  pre- 
viously. The  foil  was  taped  to  a  piece  of  graph  paper  having 
ten  divisions  per  one-half  inch.   The  foil  was  then  cut  into 
narrow  strips  approximately  0.050  inches  wide.  These  small 
strips  of  foil  were  then  placed  on  a  special  plancet  and 
counted  using  a  gas  flow  proportional  counter.  At  large  dis- 
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tances  from  the  rod,  the  small  strips  were  counted  every 
quarter  inch.   In  the  proximity  of  the  rod,  each  individual 
strip  was  counted.   At  each  position  where  a  count  was  taken, 
the  time  when  the  counting  period  began  was  recorded. 

Counting  Facilities 

The  facility  used  to  count  the  irradiated  gold  wire  was 
one  that  had  been  previously  constructed  by  the  International 
Institute  specifically  designed  to  count  irradiated  wires. 
The  apparatus  consisted  of  an  Argonne  scaler,  a  Nuclear  Chi- 
cago Model  3031B  Geiger-Mueller  counting  chamber,  and  a  wire 
guide. 

The  wire  guide  and  counting  chamber  were  so  constructed 
that  only  a  small  section  of  the  wire  could  be  "seen"  by  the 
end  window  G-M  tube.  The  portion  of  the  wire  guide  below  the 
G-M  tube  was  constructed  of  lead  to  prevent  radiation  from 
the  wire  from  being  scattered  into  the  counting  region  of  the 
tube;  see  Figure  2i|. 

The  equipment  had  been  checked  by  a  member  of  the  Labor- 
atory staff,  and  it  was  found  to  operate  best  at  a  voltage  of 
11+00  volts  with  a  discriminator  setting  of  175*  When  a  Chi- 
square  test  was  run  using  a  standard  beta  ray  source,  the 
equipment  was  found  to  be  operating  with  a  high  confidence 
probability,  QQ%. 
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To  count  the  gold  foil,  a  gas  flow  proportional  counter 
and  a  scaler  were  used.  The  scaler  was  the  Argonne  type,  and 
the  counting  chamber  was  a  lead  shielded  "pig"  type  into  which 
a  continuous  flow  of  P-10  gas,  a  mixture  of  10$  methane  and 

90%   argon,  was  passed  to  maintain  the  background  radiation 
level  as  low  as  possible. 

A  member  of  the  Laboratory  staff  had  checked  this  appar- 
atus and  found  that  it  operated  best  at  a  voltage  of  3150 
volts  with  a  discriminator  setting  of  SO*     When  a  Chi-square 
test  was  run  using  a  standard  beta  ray  source,  this  equipment 
was  found  to  be  operating  with  a  high  confidence  probability, 
Qo%. 

Treatment  of  Raw  Data 

The  data  for  each  gold  wire  and  gold  foil  were  corrected 
for  the  background  counts,  the  decay  of  the  sample  that  occurred 
previous  to  being  counted,  and  the  dead  time  of  the  counting 
equipment.  These  corrected  data  are  listed  in  Tables  2  through 

7  in  the  form  of  counts  per  minute.  The  deviation  reported 
is  the  standard  deviation  of  a  single  observation.  Each  sample 
was  counted  for  a  sufficient  length  of  time  to  reduce  the  stan- 
dard error  to  approximately  1$.  An  IBM  1620  computer  program 
was  used  to  make  the  corrections,  and  the  program  is  described 
in  Appendix  £. 
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The  data  points  near  the  outer  boundary  of  the  graphite 
stringer  were  rejected  for  the  determination  of  the  effective 
rod  radius  because  they  were  not  consistent  with  the  remainder 
of  the  data. 

Description  of  Tables 

The  tables  of  the  original  data,  Tables  2  through  7,  are 
identified  by  the  following  notation: 

Case  1.   pne  inch  rod  with  gold  wire  near  source. 

Case  2.   One  inch  rod  with  gold  wire  away  from  source. 

Case  3.  One  inch  rod  with  gold  foil  near  source. 

Case  k»     One  inch  rod  with  gold  foil  away  from  source. 

Case  5*  Two  inch  rod  with  gold  wire  near  source. 

Case  6.  Two  inch  rod  with  gold  wire  away  from  source. 

The  phrases  "near  source11  and  "away  from  source"  refer 
to  the  location  of  the  gold  foil  or  wire  with  respect  to  the 
reactor  core  during  the  period  of  irradiation.  Thus,  near 
source  means  that  the  foil  (or  wire)  was  placed  in  the  spe- 
cial graphite  stringer  in  the  position  nearest  the  core. 
Away  from  source  means  that  the  foil  (or  wire)  was  placed  in 
the  special  graphite  stringer  in  the  position  farthest  from 
the  core* 

The  data  in  these  tables  have  been  corrected  for  the 
background  counts,  the  decay  of  the  sample  that  occurred 
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previous  to  being  counted,  and  the  dead  time  of  the  counting 
equipment.  The  deviation  that  is  reported  is  the  standard 
deviation  of  a  single  observation. 
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TABLE  2. 
ORIGINAL  DATA  FOR  CASE 


DISTANCE  FROM 

COUNTS 

STANDARD 

CENTER  OF 

PER 

DEVIATION 

ONE  INCH  ROD 

MINUTE 

( INCHES) 

.00000 

10115 

102 

.06250 

10058 

101 

.12500 

10001 

101 

.18750 

10112 

102 

.25000 

10438 

103 

.31250 

10636 

104 

.37500 

10950 

106 

.43750 

11062 

106 

.50000 

11593 

109 

.56250 

11892 

110 

.62500 

12100 

111 

.68750 

12358 

112 

.75000 

12423 

112 

.81250 

12935 

115 

.87500 

12982 

115 

.93750 

13321 

116 

1.00000 

13671 

118 

1.06250 

13568 

117 

1.12500 

13658 

118 

1.18750 

13840 

118 

1.25000 

1368  2 

118 

1.31250 

13859 

118 

1.37500 

13949 

119 

1.43750 

13992 

119 

1.50000 

13958 

119 

1.56250 

14007 

119 

1.68750 

14375 

120 

1.81250 

14562 

121 

1.93750 

14762 

122 

2.06250 

14941 

123 

2.18750 

14768 

122 

2.31250 

14894 

122 

2.43750 

15080 

123 

2.56250 

15188 

124 

2.68750 

14766 

122 

2.81250 

14617 

121 

2.93750 

14260 

120 

3.06250 

14167 

119 

3.18750 

14060 

119 

3.31250 

14086 

119 

3.43750 

13980 

118 
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TABLE  3. 
DATA  FOR 


CASE  2 


DISTANCE  FROM 

COUNTS 

STANDARD 

CENTER  OF 

PER 

DEVIATION 

ONE  INCH  ROD 

MINUTE 

( INCHES) 

.03125 

8434 

94 

.09375 

8387 

94 

.15625 

8613 

95 

.21875 

8783 

96 

.28125 

8921 

97 

.34375 

9346 

99 

.40625 

9744 

101 

.46875 

10325 

104 

.53125 

10585 

105 

.59375 

11045 

107 

.65625 

11114 

108 

.71875 

11754 

110 

.78125 

11938 

111 

.84375 

12039 

112 

.90625 

12158 

112 

.96875 

12164 

112 

1.03125 

12326 

113 

1.15625 

12740 

115 

1.28125 

12775 

115 

1.40625 

12807 

115 

1.53125 

12916 

115 

1.65625 

12982 

115 

1.78125 

13033 

116 

1.84375 

13406 

117 

1.90625 

13234 

116 

1.96875 

13004 

115 

2.03125 

13189 

116 

2.15625 

13188 

116 

2.28125 

13508 

117 

2.40625 

13545 

117 

2.53125 

13863 

119 

2.65625 

13660 

118 

2.78125 

13257 

116 

2.90625 

13230 

116 

3.03125 

13174 

116 

3.15625 

13179 

116 

3.28125 

13161 

115 
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TABLE  4. 
ORIGINAL  DATA  FOR  CASE  3 


DISTANCE  FROM 

COUNTS 

STANDARD 

CENTER  OF 

PER 

DEVIATION 

ONE  INCH  ROD 

MINUTE 

( INCHES) 

.02500 

9677 

46 

.07500 

9714 

48 

.12500 

9739 

49 

^17500 

10016 

48 

.22500 

10092 

47 

.27500 

10293 

47 

.32500 

10996 

59 

.37500 

11552 

49 

.42500 

11947 

56 

.47500 

12470 

57 

.52500 

13268 

63 

.57500 

13660 

60 

.72500 

14317 

64 

.87500 

14796 

67 

1.02500 

15031 

62 

1.17500 

15177 

62 

1.32500 

15374 

70 

1.47500 

15737 

70 

1.62500 

15797 

66 

1.77500 

16105 

75 

1.92500 

16307 

66 

2.07500 

16260 

69 

2.22500 

16216 

70 

2.37500 

16024 

64 

2.52500 

16604 

78 

2.67500 

16145 

69 

2.82500 

16651 

73 

2.97500 

16622 

67 

3.12500 

16636 

71 

3.27500 

16398 

64 

3.42500 

16706 

68 

3.57500 

16916 

77 

3.77500 

16315 

63 

3.92500 

16356 

71 
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TABLE  5. 
DATA  FOR 


CASE  4 


DISTANCE  FROM 

COUNTS 

STANDARD 

CENTER  OF 

PER 

DEVIATION 

ONE  INCH  ROD 

MINUTE 

(  INCHES) 

.02500 

5114 

29 

•07500 

5019 

20 

.12500 

4857 

25 

.17500 

4852 

20 

.22500 

5205 

30 

.27500 

5024 

22 

.32500 

5025 

22 

.37500 

5157 

23 

.42500 

5194 

26 

.47500 

5287 

26 

.52500 

5441 

24 

.62500 

5766 

29 

.72500 

5780 

25 

.82500 

5975 

26 

.92500 

6079 

31 

1.02500 

5988 

26 

1.17500 

6293 

28 

1.32500 

6276 

28 

1.47500 

6248 

31 

1.67500 

6340 

32 

1.82500 

6396 

29 

1.97500 

6875 

30 

2.12500 

6020 

30 

2.27500 

6743 

31 

2.42500 

6781 

39 

2.57500 

6789 

34 

2.72500 

6682 

30 

2.87500 

6710 

33 

3.02500 

6598 

30 

3.17500 

622  5 

31 

3.32500 

6443 

36 

3.47500 

6707 

31 

3.62500 

6825 

31 

3.77500 

6165 

29 

3.92500 

6324 

26 
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TABLE  6. 
ORIGINAL  DATA  FOR  CASE  5 


DISTANCE  FROM 

COUNTS 

STANDARD 

CENTER  OF 

PER 

DEVIATION 

TWO  INCH  ROD 

MINUTE 

(  INCHES) 

.03125 

4555 

48 

.09375 

4805 

50 

.15625 

4867 

50 

.21875 

5079 

51 

.28125 

5330 

52 

.34375 

5442 

53 

.40625 

5754 

54 

.46875 

6022 

56 

.53125 

6163 

56 

.59375 

6427 

57 

.65625 

6779 

59 

.71875 

7273 

71 

.78125 

7640 

72 

.84375 

8035 

74 

.90625 

8508 

76 

.96875 

8589 

77 

1.03125 

8962 

78 

1.09375 

9019 

79 

1.15625 

9026 

79 

1.21875 

9225 

80 

1.28125 

9354 

80 

1.34375 

9223 

80 

1.40625 

9512 

81 

1.46875 

9621 

81 

1.53125 

9745 

82 

1.59375 

10079 

102 

1.65625 

10580 

105 

1.71875 

10556 

105 

1.78125 

10938 

106 

1.84375 

11198 

108 

1.90625 

11360 

109 

1.96875 

11631 

110 

2.09375 

11719 

110 

2.21875 

11569 

110 

2.34375 

11299 

108 

2.46875 

11202 

108 

2.59375 

11271 

108 

2.71875 

11392 

109 

2.84375 

11782 

111 

2.96875 

11820 

111 

3.09375 

12374 

113 

3.21875 

12529 

114 
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TABLE  6.  CONTINUED 
ORIGINAL  DATA  FOR  CASE  5 


DISTANCE  FROM 

COUNTS 

STANDARD 

CENTER  OF 

PER 

DEVIATION 

TWO  INCH  ROD 

MINUTE 

( INCHES) 

3.34375 

12462 

114 

3.46875 

12799 

115 

3.59375 

13211 

117 

3,71875 

13599 

119 

3.84375 

13865 

120 

3.96875 

13446 

118 

4.09375 

13554 

119 

4.21875 

13132 

117 

4.34375 

12289 

113 

4.46875 

12119 

113 

4.59375 

12134 

113 

4.71875 

11941 

112 

4.84375 

12627 

115 

4.96875 

12103 

113 

5.09375 

11778 

111 

5.21875 

11621 

110 

5.34375 

12122 

113 

5.46875 

11979 

112 

5.59375 

13159 

118 

5.71875 

12965 

117 

5.84375 

12845 

116 

5.96875 

12249 

113 

6.09375 

11762 

111 

6.21875 

'  11543 

110 
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TABLE  7. 
ORIGINAL  DATA  FOR  CASE  6 


DISTANCE  FROM 

COUNTS 

STANDARD 

CENTER  OF 

PER 

DEVIATION 

TWO  INCH  ROD 

MINUTE 

( INCHES) 

.OOOOO 

3232 

43 

•06250 

3301 

43 

.12500 

3303 

43 

.18750 

3465 

44 

.25000 

3613 

45 

.31250 

3729 

46 

437500 

3957 

47 

•43750 

4134 

48 

.50000 

442  7 

49 

.56250 

4638 

50 

.62500 

4792 

51 

.68750 

5074 

53 

.75000 

5219 

62 

.81250 

5620 

64 

.87500 

5963 

65 

.93750 

6424 

74 

1.00000 

6806 

76 

1.06250 

7147 

78 

1.12500 

7141 

78 

1.18750 

7582 

80 

1.25000 

7644 

81 

1.31250 

7873 

82 

1.37500 

8135 

83 

1.43750 

8029 

82 

1.50000 

8164 

83 

1.56250 

8112 

93 

1.62500 

8208 

93 

1.68750 

8348 

94 

1.75000 

8408 

94 

1.81250 

,8401 

94 

1.87500 

8432 

94 

1.93750 

8480 

94 

2.00000 

8510 

94 

2.06250 

8784 

96 

2.18750 

9154 

98 

2.31250 

9048 

97 

2.43750 

9279 

98 

2.56250 

9392 

99 

2.68750 

9229 

98 

2.81250 

9617 

100 

2.93750 

9734 

100 

3.06250 

9644 

100 
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TABLE  7.  CONTINUED 
ORIGINAL  DATA  FOR  CASE  6 


DISTANCE  FROM 

COUNTS 

STANDARD 

CENTER  OF 

PER 

DEVIATION 

TWO  INCH  ROD 

MINUTE 

(  INCHES) 

3.18750 

10159 

103 

3.31250 

10230 

103 

3,43750 

10158 

102 

3.56250 

10593 

105 

3.68750 

10667 

105 

3.81250 

10848 

106 

3.93750 

10693 

105 

4.06250 

10443 

104 

4.18750 

10195 

102 

4.31250 

10518 

104 

4.43750 

11081 

107 

4.56250 

11202 

107 

4.68750 

11697 

109 

4.81250 

11512 

108 

4.93750 

11884 

110 

5.06250 

11852 

110 

5.18750 

11562 

109 

5.31250 

11152 

107 

5.43750 

10988 

106 

5.56250 

10576 

104 

5.81250 

10664 

104 

6.06250 

10477 

103 

6.18750 

10443 

103 

6.31250 

10254 

102 

6.43750 

10185 

102 

6.56250 

9948 

100 

6.68750 

9791 

100 
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APPENDIX  E 


IBM  1620  Program  to 
Correct  Raw  Data 


This  computer  program  was  written  to  correct  the  raw 
counted  data  for  background  counts,  the  counter  dead  time, 

and  the  decay  of  the  sample  previous  to  being  counted.  The 
program  was  written  in  FORTRAN  I,  and  it  utilized  the  auto- 
matic floating  point  subroutines  EXP  and  SQRT.  The  logic 
diagram  and  the  object  program  are  given  in  this  section. 

Input  for  this  program  consisted  of  the  experimental 
value  of  the  counts,  the  length  of  time,  in  minutes,  that 
elapsed  to  record  the  counts,  the  value  of  the  exponent  for 
the  decay  correction,  and  the  number  of  the  sample  being 
corrected.  A  base  value  was  also  included  to  permit  the 
counts  to  be  based  on  any  desired  length  of  time.  The  data 
for  each  point  was  punched  on  a  single  card  according  to  the 
specified  format  statement  and  read  into  the  computer  after 
the  object  program  and  the  subroutines  had  been  stored. 

The  output  can  be  in  either  punched  or  typed  form  depend- 
ing on  the  setting  of  sense  switch  two.   One  output  card  was 
punched  for  each  input  card  that  was  read.  The  output  con- 
sisted of  the  corrected  count  rate,  the  number  of  the  sample 
that  was  corrected,  and  the  standard  deviation  of  the  count 
rate. 
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The  operating  time  for  this  program  was  approximately 
one  second  to  correct  one  data  point. 

To  correct  the  raw  data  counts  recorded  using  the  gold 
foil,  a  separate  program  was  written.   This  program  accounted 
for  the  variation  in  the  weights  of  the  individual  foils. 

The  input  for  this  program  was  identical  to  that  of  the 
program  that  processed  the  wire  data  except  the  weight  of 
the  foil,  in  milligrams,  was  included.   The  data  for  each 
point  are  punched  on  a  single  card  according  to  the  specified 
format  statement  and  read  into  the  computer  after  the  object 
program  and  the  subroutines  had  been  stored. 

The  output  for  this  program  was  identical  to  that  of  the 
wire  data  processing  program.  The  operating  time  for  this 
program  was  approximately  one  second  per  data  point. 


Table  8.   Input  Data  For  the  IBM  1620  Raw  Wire 
Data  Correction  Program 


Symbol  Explanation 


BASE      Length  of  time  on  which  the  count  rate  is  based 

in  minutes 
CR        Total  number  of  counts  recorded  for  the  sample 
M        Length  of  the  counting  period  in  minutes 
E        Exponent  for  the  decay  correction,  At 
N         Sample  Number 
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Table  9.   Input  Data  For  the  IBM  1620  Raw  Foil 
Data  Correction  Program 


Symbol  Explanation 


CR  Total  number  of  counts  recorded  for  the  sample 

M  Length  of  the  counting  period  in  minutes 

E  Exponent  for  decay  correction,  At 

N"  Sample  number 

WT  Sample  weight  in  milligrams 
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0  PROGRAM   FOR  PROCESSING   THE  WIRE  DATA 

1  FORMAT    (E10.4,I5,E10.4.I5) 

2  FORMAT   (E14.7,I6,E14.7) 
B=28. 3833333 

READ1,    BASE 

20  READ1,CR,M,E,N 
TIME=M 
CPM=CR/TIME 

CCRPM=CPM/(1.-CPM*4. E-6 ) 
CON=BASE*EXP(E) 
CNT=(CCRPM-B)*CON 

DEV=SQRT((CCRPM+B)*TIME)*CON/TIME 
IF( SENSE  SWITCH  2)21,22 

21  PRINT2,CNT,N,DEV 

22  PUNCH2,CNT,N,DEV 
GO    TO   20 

END 


0    PROGRAM  FOR  PROCESSING  RIBBON  DATA 

1  FORMAT   (E10.4,I5,E10.4.I5,E10.4) 

2  FORMAT   (E14.7,I6,E14.7) 
SIGWS=6.25E-4 

B=15.0 

20  READ1,CR,M,E,N,WT 
TIME=M 
CPM=CR/TIME 

ENT=CPM/(l.-CPM*6.745E-8) 
EN=(ENT-B)*TIME 

CNT=2. *EN*EXP(  E)/(WT*TIME) 

SIGNS=(ENT+B)*TIME 

DEV=CNT*SQRT( SIGNS/( EN**2 )+SIGWS/(WT**2 ) ) 

IF( SENSE  SWITCH  2)21,22 

21  PRINT2,CNT,N,DEV 

22  PUNCH2,CNT,N,DEV 
GO   TO  20 

END 


99 


APPENDIX  P 

Description  and  Explanation  of  the  IBM 
1620  Computer  Program  to  Determine  the 
Effective  Rod  Radius  Prom  a  Set  of 
Experimental  Data 

To  determine  the  effective  rod  radius  from  the  experi- 
mental data,  it  is  necessary  to  select  a  method  to  calculate 
the  coefficients  in  the  flux  equations 

<tf  (r)   =T  A^  XJ/3^  12)  (P-i ) 

and 

where  3  is  summed  over  the  number  of  harmonics  included  in 
the  analysis  of  the  data,  and  i  is  the  position  at  which  the 
flux  was  measured.   To  calculate  these  coefficients,  the 
method  of  weighted  least  squares  is  used. 

A  weightlngt  factor  in  a  least  squares  problem  may  be 
utilized  in  two  ways:   one,  by  requiring  that  the  sum  of  the 
weighted  residuals  squared  be  a  minimum;  or  two,  by  requiring 
that  the  residual  between  the  calculated  value  at  any  point 
should  lie  within  the  statistical  standard  deviation  of  the 
experimental  measurement.  A  weighting  factor  of  unity  satis- 
fies the  first  method  best.   Scarborough  (14)  points  out  that 
the  residuals  squared  are  best  weighted  with  the  reciprocal 
of  the  standard  deviation,  1/0~2  ,  if  the  second  method  is  to 
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be  satisfied. 

Recall  from  the  theory  section  that  the  constants  or- 
and  yS-  are  related  by  the  equations 

cr 

/3j  -  y; ->£-*£  (F-3) 


where  T=  J0(.ar-Rl)/\i(^-(?l)and  >C  and  K^"  are  constants 
of  the  rod  and  the  moderator. 

Thus,  it  is  evident  that  for  a  given  value  of  Rl  and 
R2,  the  simultaneous  solution  of  Eqs.  (F-3)  and  (F-l|)  will 
yield  a  proper  value  for  <%.■    and  A/  .   In  order  to  begin  the 
calculation  with  a  reasonable  value  for  oc-    and  Aj   ,  an 
initial  value  for  o(j    is  estimated  by  using  the  boundary  con- 
dition  that  the  flux  goes  to  zero  at  the  outer  boundary,  R2. 

Once  the  value  of  oCj  and  Aj    for  a  given  Rl  is  deter- 
mined, the  unknown  coefficients  in  the  flux  equations  are 
calculated  by  the  method  of  least  squares.   The  equation  that 
must  be  solved  is 

E  =YWii<j*  -E/V2  (F-5) 

where 

M^   =  A^IoC/SjTt)  0<r«R1  (P-6) 

M^     =  &j  J.C^-yv)  +  CjY.(«jn)  R1<rx<  Rl  ,  (F-7) 
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In  order  that  this  method  might  fully  be  utilized,  it  is 
necessary  to  express  all  the  unknown  coefficients  in  terms 
of  one  of  their  number.   The  ooeffioient  Bj  is  seleoted  as 
the  main  unknown.   Using  the  boundary  condition  that 


<£M(R2)  -0 
the  expressions  for  0.  and  A.  are  obtained: 

where  T^  =  JQ(  <*£  R2)/YQ(  o£  R2). 

Operating  on  Eq..  (F-5)  in  the  usual  manner,  the 
equation  for  the  determination  of  the  unknown  coeffi- 
cients is  obtained: 


(F-8) 
(F-9) 


(F-10) 
(F-ll) 


£BiI^M^MKjL  -rv^M^ 


(F-12) 


where 


M^^UKR^-T.xc^mn/x^^)}1^^   °^RA  (p"13) 


HKi=    J.(^k^)-TkY0Kiq:) 


RA<  r«R2 


and  as  before 


(F-14) 


TK   -    J0C^R2)/Yo(^KR2). 


(F-15) 
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After  the  coefficients  A.,  B,,  and  Cj  are  obtained,  the 
least  squared  error  is  calculated  using  Eq.  (P-5) .   This  value 
is  stored,  Rl  is  incremented,  and  the  procedure  is  repeated, 

starting  with  the  determination  of  ci*    and  Aj    . 

The  value  of  Rl  is  incremented  twice  from  the  original 
value,  and  the  least  squared  error  at  each  of  these  values 
of  Rl  is  subtracted  in  the  manner  of  forward  differences. 
These  differences,  along  with  the  average  value  of  the  radii 

used  to  obtain  the  errors  are  then  used  to  determine  a  straight 
line  that  intersects  the  axis  of  the  graph,  A  Error  vs.  Rl, 
at  the  position  of  an  estimate  of  the  minimum  on  the  error 
curve.   (See  Figure  25) 

Using  this  value  of  Rl  as  the  best  estimate  of  the  mini- 
mum of  the  error  curve,  the  least  squared  error  for  this 
position  is  calculated  and  stored,  Rl  is  incremented,  once 
in  the  direction  of  increasing  Rl  and  once  in  the  direction 
of  decreasing  Rl,  and  the  least  squared  error  at  each  of  these 
positions  is  calculated,  stored,  and  used  to  determine  a 
A  Error  and  a  new  Rl  as  before. 

After  the  error  is  calculated  with  this  second  estimate 
of  the  minimum  of  the  error  curve,  it  is  compared  to  the  error 
calculated  with  the  first  estimate  of  the  minimum.   If  these 
errors  are  close  together  and  the  absolute  value  of  their 
difference  is  less  than  or  equal  to  a  given  criterion,  the 
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calculation  Is  complete.   Then  the  values  of  the  unknown 

coefficients,  o(:    and  A-    ,  the  value  of  the  effective  rod 
radius,  and  the  flux  plot  that  results  In  the  minimum  least 
squared  error  is  punched. 

If  the  error  criterion  Is  not  satisfied,  the  two  radii 
that  estimate  the  minimum  error  are  compared.   If  they  are 
close  together  and  the  absolute  value  of  their  difference  is 
less  than  or  equal  to  a  preset  value,  the  calculation  is 
complete,  and  the  required  Information  is  punched.   If  the 
radii  accuracy  criterion  is  not  satisfied,  the  errors  associ- 
ated with  these  radii  are  compared  to  determine  which  is  the 
smaller.  The  smallest  error  and  the  associated  radius  is 
stored,  and  the  iteration  to  determine  the  minimum  error  pro- 
ceeds. 

The  amount  of  output  obtained  during  the  iterations,  and 
the  type  of  analysis  to  be  employed  during  the  calculations 

can  be  controlled  by  the  proper  setting  of  the  four  sense 
switches.   If  sense  switch  1  is  ON  while  the  input  data  is 
being  read  into  the  computer,  the  initial  values  of  the  count 
rate,  the  position,  with  respect  to  the  center  of  the  rod,  at 
which  the  count  was  recorded,  and  the  weighting  factors  are 
punched.   If  switch  1  is  OFF,  all  of  this  output  is  omitted. 
When  switch  1  is  ON  while  the  iterations  are  in  progress, 
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the  intermediate  values  of  c*v  and  /#'  from  the  solution  of 
Eq.  (F-4)  are  punched.   This  output  is  omitted  if  switch  1 
is  turned  OFF  during  the  iterations. 

When  sense  switch  2  is  turned  ON  during  the  iterations, 
the  value  of  the  least  squared  error  and  the  radius  at  which 
it  was  calculated  is  punched.  If  switch  2  is  OFF,  this  out- 
put is  eliminated.  Turning  switch  2  ON  as  the  input  data  is 
being  read  causes  the  transcendental  equation  for  X  in  the 
rod  to  be  solved.  If  the  switch  is  OPP  as  the  data  is  read, 
the  diffusion  theory  value  for  K  is  calculated  and  used  for 
the  determination  of  the  effective  rod  radius. 

Turning  sense  switch  3  ON  will  cause  an  exit  from  the 
program  during  the  iterative  calculation  for  o(:    and  /3y  or 
following  the  punching  of  the  least  squared  error  and  Rl. 
Sense  switch  2  must  be  turned  ON  to  exit  following  the  punch- 
ing of  the  error  and  Rl.  When  this  exit  occurs,  all  of  the 
output  that  results  when  the  calculation  ends  in  a  normal 
fashion  is  punched.   Oare  must  be  exercised  when  this  exit 
is  chosen.  A  minimum  of  eight  iterations  must  have  occurred 
before  this  type  of  exit  is  attempted.   If  sense  switch  3  is 
OPP,  the  exit  is  omitted,  and  the  end  of  calculation  type 
output  is  eliminated. 

The  internal  subroutine  CRAM  also  operates  on  sense 
switch  3.   If  switch  3  is  ON  and  the  calculation  of  the 
unknown  flux  coefficients  is  in  process,  the  value  of  the 
matrix  elements  (Eq.  (P-12))  is  typed  before  the  matrix  is 
solved.   When  switch  3  is  OPP,  this  output  is  omitted. 
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Because  of  the  nature  of  the  operations  connected  with  the 
use  of  sense  switch  3,  It  must  be  used  with  caution. 

If  the  program  is  compiled  with  the  traoe  Instructions 
included,  sense  switch  4  controls  typed  output.   The  compi- 
lation of  the  trace  instructions  and  the  use  of  sense  switch 
4  are  primarily  for  the  "debugging"  of  new  programs.   When 
switch  4  is  ON,  the  values  to  the  left  of  all  equal  signs 
are  typed  out.   This  type  of  output  greatly  slows  the  oper- 
ation of  the  computer.   So,  unless  a  program  is  troublesome, 
sense  switch  4  is  normally  left  turned  OFF. 
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Table  10.   Input  Data  for  the  IBM  1620  Effective 
Rod  Radius  Program 


Symbol 


Explanation 


SIGR      Total  macroscopic  cross  section  of  the  rod,  cm" 
SIGSR     Macroscopic  scattering  cross  section  of  the  rod,  cm 
SIGAR     Macroscopic  absorption  cross  section  of  the  rod,  cm 
SIGAG     Macroscopic  absorption  cross  section  of  the 

moderator,  cm"1 
SIGSG     Macroscopic  scattering  cross  section  of  the 

moderator,  cm""" 
AMTJZ      Average  cosine  of  the  scattering  angle  in  the  rod 
AMTJ1      Average  cosine  of  the  scattering  angle  in  the 

moderator 
Rl        Radius  of  the  rod,  inches 
R2        Outer  boundary  radius  of  the  moderator,  inches 

N2  Number  of  data  points 

Jl  Number  of  harmonics  used  in  the  analysis 

IVD  Divisor  to  convert  Kl  to  a  distance  from  center 

of  rod,  inches 

Kl  Location  of  sample  from  center  of  rod 

PHI  Corrected  count  rate  of  sample,  cpm 

M  i           Sample  number 

WT  Standard  deviation  of  the  corrected  count  rate,  cpm 


-1 

-1 


108 


o 

o 

81 

Calculate 
Initial 

a> 

°H- 

o   *i 

o 

o 

J 

, 

a> 

•■     *i 

O    Q 

3 

i>  o" 

o«  J 

o  x 

CO     -    - 

„cr  o  — 

^°cr  co  co  => 

£0002 

i  cr  co  co  co  < 


o- 
ii  ii 

-IS 


a> 
CO 


^^ 

^^ 

<fr 

2 

cr 

K 

o 

ERRO 
=  Rl 

a: 
cr 

UJ 

ii 

ii  •-» 

^^ 

-» x: 

*z 

•— i «— < 

5q 

a: 

_o 

or-i 

ii  * 

CT  O 

a: 

uj  <r 

2iu 

<T> 

^ 

•  -cr> 

•J"6' 

a) 

fc  u. 

> 

«-    & 

o 

O    m 

CO 

Lu  ZZ 

C  "3 

feetJ 

£uj  h- 


CVJ 

(0 

CJ  o 

II  " 

♦»         _ 

(DmO 

COl^O 

Kw 


cr 


".  ..  cr_9.  "  2 


"S  <0  "  »-<  o 
co  cr  cr2crcro 


3?tj     o 

_  o  cr  <i» 
02  ocr 
co~ 


ro 

a>  _       — 

*■•             ^™ 

0  cr  +  0 

3         2»" 

0  *  ,. 

o«"o 

0220 

il 

X 

*-             c 

0  -a        0 

2  a>  H  — 

£2" 

«T3    O-O 

>  0  cr  « 

02  0  cr 

G0~ 

1 

X 

k. 

O     , 

2i2~ 

C   "3 

■»-  a>     " 

!S  E  ~ 

Cuh 

■D 

O) 


O 


ico  w 


CO 

3 

O 

cr 

"O 

o 

cr 

a> 

> 

.  — 

o 

Q) 

«*- 
•♦- 

UJ 

o 

c 


0) 
CD 

Q 


£ 
o 

o 


o 

OJ 

CD 


en 


E 
o 

o 


o 


109 


O 

<z 

•o 
o 
q: 

o 

> 

o 

CD 


<L> 

C 

£ 

_ 
o 

o 
Q 


o 
o 


O 
(0 


CD 


< 

1j     crfo 

O         li 

±  +  s£ 

/2Sh 

0>J1  "  <  o 

cocr  2  cro 

£ 

a 

u. 
en 
a 


o 


110 


C      EFFECTIVE  ROD  RADIUS  PROGRAM      WW  PCRATH       5/1/63 

DIMENSION  TA( 5,6) »B<  5) ,PHI  (70)  »RA( 70) ,WW( 70) 

DIMENSION  AIZ(5,30)  ,AJZ(5»70) ,AYZ( 5.70 )  »P(  5*70) 

DIMENSION  W(5)»XX(5),YY(5)»C(5),A(5) ,R0LD1(5) 

DIMENSION  AA( 5)  tBG( 5 ) . ALPHA ( 5)»T(5)»A1A(5) 

DIMENSION  RAD(5) , ERROR (5) ,ERR(5) »DERR(5) ,CONST(5) 

G  =  0. 

N  =  l 

I1  =  G 

READ5»SIGR,SIGSRtSIGAR,SIGAG»SIGSG 

READ5.  AMUZ,  AMU1 

READ5»CRIT»DELTA.VAR 
C      INPUT  ABOVE  IS  3  CARDS 

DZER=l./( 3.*(SIGAR+SIGSR*( l.-AMUZ) ) ) 

DONE=l./( 3.*(SIGAG+SIGSG*( 1.-AMU1) ) ) 

AKP1=(SIGAG/DGNE)*(2.54)**2 

AKP09=(SIGAR/DZER)*(2.54)**2 

AAKPO=SQRT(SIGAR/DZER) 

AKP01=AAKP0+l.E-4 
C      THIS  SENSE  SWITCH  SELECTS  THE  TYPE  OF 
C      ANALYSIS  FOR  KAPPA  OF  ROD 

IF(SENSESWITCH2)29,2 
2  DZER=1./(3.*SIGSR) 

DONE=l./( 3.*SIGSG) 

AKP0=(SIGAR/DZER)*(2.54)**2 

AKP1=(SIGAG/D0NE)*(2.54)**2 

GO  TO  15 
C      TRANSCENTENTAL  SOLUTION  FOR  KAPPA  IN  ROD  TO  STATEMENT  15 

29  TEST2=AAKP0-AKP01 
IF(TEST2)28»33»26 

28  TEST2=-TEST2 

26  IF(TEST2/AAKPO-CRIT)33»33»27 

27  SOME=EXP( AAKPO/SIGSR) 
S0ME=AAKPO/SIGR-(SOME-(1.0/SOME) )/(SOME+( 1.0/SOME) ) 
IF(S0ME)36. 33.36 

36  IF( 11)30,30,31 

30  I1=N 
AKP01=AAKP0 
AAKPO=AAKPO+DELTA 
S0ME1=S0ME 

GO  TO  29 

31  IF(AAKPO-AKP01)32»33,32 
33  AKP0=(AAKP0*2.54)**2 

GO  TO  15 

32  SLOP= ( S0ME-S0ME1 ) / ( AAKP0-AKP01 ) 
IF(SL0P)37,38,37 

38  PUNCH95 
GO  TO  30 

37  AKP01=AAKP0 
S0ME1=S0ME 

AAKPO=AAKPO-SOME/SLOP 
GO  TO  29 


Ill 


15  READ5,  Rl»  R2»  ACC »  DL»  DLTA»VAL 
READ6.N2»J1,IVD 
C      INPUT  ABOVE  IS  2  CARDS 
DIV=IVD 
DC  35  1=1, N2 
READ6.K1 
AB  =  K1 

RA(  I  )=AB/DIV 
IF(RA( I )-Rl) 34,34,35 

34  NRCD=I 
L1=I+1 

35  CONTINUE 

C      INPUT  ABOVE  IS  N2  CARDS 

DO  41  1  =  1, N2 

READ  99,  PHI { I ) ,M,WT 
41  WW(  I  )=1./(WT*WT) 
C      INPUT  ABOVE  IS  N2  CARDS 

D2ER=DZER/2.54 

D0NE=D0NE/2.54 
C      THE  FOLLOWING  7  CARDS  PUNCH  INPUT  DATA 

PUNCH 1, ACC, AKPO,AKPl»DIV 

PUNCH 1,DZER, DONE, DL,AKP09 

PUNCH 1,R1,R2,DLTA,VAL 

PUNCH 1,CR IT, DELTA, VAR,SIGR 

PUNCH 1, SI  GSR, SI  GAR, SI GAG, S 1 6SG 

PUNCH1,AMU2,AMU1 

PUNCH6,N2,J1 
C      THIS  SENSE  SWITCH  IS  FOR  THE  PUNCHING  OF  R,PHI,WT 

IF(SENSE  SWITCH  1)504,506 
504  PUNCH994 

DO  980  I=l,NROD 

980  PUNCH1,RA( I ) ,PHI ( I ) ,WW( I ) 
PUNCH993 

DO  981  I=L1,N2 

981  PUNCH1,RA( I ) ,PHI( I ) ,WW( I ) 
506  ALPHAd  )=2. 4048256 

ALPHA(2)=5. 5200781 

ALPHA(3)=8. 6537279 

ALPHA(4)=11. 7915344 

ALPHA(5 )=14. 9309177 

L  =  G 

M  =  N 

KI=N 
C      THIS  LOOP  CALCULATES  THE  INITIAL  VALUES 
C      OF  ALPHA  AND  BETA 

DO  17  1=1, Jl 

AA( I )=ALPHA( I )/R2 

A1A( I )=AA( I )+l.E-4 

ARG1=AKPC-AKP1-AA( I ) **2 
17  BG( I )=SQRT(ARG1 ) 

I  I  I  =G 
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C      THIS  LOOP  SOLVES  THE  TRANSCENDENTAL  EQUATION 
C      FOR  ALPHA  AND  BETA 

13  DO  16  1  =  1. Jl 
KKK  =  G 

C      THIS  SENSE  SWITCH  IS  TC  EXIT  FROM  THE  PROGRAM  BEFORE 
C      THE  CALCULATION  IS  COMPLETE 

10  IF(SENSESWITCH3)71 .11 

11  W< I )=BG( I )*R1 
XX(  I  )=AA(  I )*R1 
YY( I )=AA( I )*R2 

T( I  )=JZER(YY(  I))/YZER(YY(I)) 

CONST ( I )= ( JZER(XX( I ) )-T( I )*YZER(XX( I ) ) )/IZER(W( I ) ) 

ANUM=DONE*AA( I )*( JONE(XX( I ) )-T( I )*YONE(XX( I ) ) ) 

ADEN=DZER*IONE(W( I ) ) *CONST ( I ) 

SUM=BG( I )+ANUM/ADEN 

TEST1=( A1A( I )-AA( I ) )/AA(I ) 

IF(TESTl) 12.316*14 

12  TEST1=-TEST1 

14  IF(TEST1-ACC)316.316.298 

298  IF(SUM)299.316.299 

299  IF( I  I  I )300, 300.304 

300  I  I  I =N 

A1A( I )=AA( I ) 
AA( I )=AA( I )+DL 

254  SUM11=SUM 
ARG1=AKP0-AKP1-AA( I ) **2 
BG( I )=SQRT( ARG1 ) 

GO  TO  10 
304  IF(AA( I )-AlA( I ) )250»316.250 
2  50  SL0PE=(SUM-SUM11)/(AA( I )-AlA< I  >  ) 

IF(SL0PE)261»260»261 

260  PUNCH96 
GO  TO  316 

C      THIS  SENSE  SWITCH  IS  FOR  PUNCHING  INTERMEDIATE 
C      VALUES  OF  ALPHA  AND  BETA 

261  IF(SENSESWITCH1)253.252 
253  IF(KKK)232. 232.223 

232  KKK=N 

PUNCH117 
223  PUNCH1.SUM. SLOPE. SG(  I )  .AA(  I ) 
252  A1A(  I  )=AA(  I  ) 

AA( I )=AA( I )-SUM/SLOPE 

GO  TO  254 
316  IF( 1-1 )255,255.16 

255  PUNCH  117 

16  PUNCH1.SUM»SL0PE.3G( I ) .AA( I ) 
C      THIS  LOOP  CHECKS  ON  THE  POINTS  IN  THE  ROD  AND  MODERATOR 

DO  21  I=1.N2 

IF(RA( I )-Rl)22. 22.21 
22  NROD=I 

L1=I+1 
21  CONTINUE 
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J2=J1+1 
C      THIS  LCCP  ZEROS  THE  MATRIX  TA<I,J) 
DC  190  J=l,5 
DC  190  1=1,6 

190  TA( J,I )=G 

IF( Jl-5)189,20l,15 
C      THIS  LCCP  PUTS  CNES  IN  TA(Itl)  AND  TA(I,6) 
189  DC  200  I=J2,5 
TA( I ,6)=N 

200  TA( I.I )=N 

C      THIS  LCCP  CALCULATES  THE  MATRIX  ELEMENTS  TA(I,J) 

201  DC  194  K=1.J1 
DC  194  J=1,K 
DC  194  1  =  1, N2 

IF< I-NRCD) 191,191,192 

191  IF( J-l)208,207,208 

207  AIZ(K» I )=IZER(BG(K)*RA( I ) ) 
TERM1=WW( I )*CCNST(K)*AIZ(K,I )*PHI ( I ) 

208  TERM2=WW( I ) *CCNST ( K ) *A I Z ( K, I )*CCNST ( J ) *AI Z ( J ♦ I ) 
GC  TC  193 

192  IF(J-1)209, 206,209 

206  AJZ(K,I )=JZER(AA(K)*RA< I ) ) 
AYZ(K»I )=YZER(AA(K)*RA( I ) ) 
P(K,I )=AJZ(K,I )-T(K)*AYZ(K,I ) 
TERM1=WW( I )*P(K»I )*PHI ( I) 

209  TERM2=WW( I )*P(K,I )*P(J,I ) 

193  IF( J-l)216,215»216 

215  TA(K,6)=TA(K,6)+TERM1 

216  TA(K,J)=TA(K,J)+TERM2 

194  TA(J,K)=TA(K,J) 

C      THIS  SENSE  SWITCH  PUNCHES  THE  VALUES  CF 
C      THE  MATRIX  ELEMENTS  TA(ItJ) 
IF(SENSESWITCH4)241,2  40 

240  DC  242  1  =  1  ,J1 

242  PUNCH  1,TA( I  ,1) ,TA( 1,2) »TA( I ,3) »TA( 1,6) 

241  IF( Jl-1 )211,211,212 

211  B(l  )=TA(1,6)/TA(1»1) 
GC  TC  213 

C      THIS  STATEMENT  SCLVES  THE  MATRIX  USING  THE  CRAM  SUBROUTINE 

212  B(l  )=CRAM(5. ) 

213  DC  74  1=1, Jl 

C( I  )=-B(  I  )*T( I ) 
A( I )=B( I )*CCNST ( I ) 
IF(SENSESWITCH4)74,2  05 
205  PUNCH1,A( I ) ,B( I ) ,C( I ) 
74  CONTINUE 
C      THIS  LOOP  CALCULATES  THE  LEAST  SQUARED  ERROR 
ERROR(M)=G 
DC  196  I=1,N2 
SUMM=G 
DC  195  J=1,J1 
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IF( I-NRCD) 197,197.198 

197  S  =  A( J)*AIZ( J,  I  ) 
GO  TO  195 

198  S  =  B( J)*AJZ( J,  I  )+C( J)*AYZ( J,  I  ) 

195  SUMM=SUMM+S 

196  ERRCR(M)=ERRCR(M)+WW( I )*(PHI ( I )-SUMM)**2 

C      THIS  SENSE  SWITCH  IS  FOR  THE  PUNCHING  OF  THE  ERROR  AND  Rl 

IF(SENSESWITCH2) 120,119 
120  PUNCH100,ERRCR(M) ,R1 

PUNCH6,M 
C      THIS  SENSE  SWITCH  IS  TO  EXIT  FROM  THE  PROGRAM  BEFORE 
C      THE  CALCULATION  IS  COMPLETE 

IF(SENSESWITCH3)71,119 
119  IF(M-4)122, 131,122 
C      THIS  SECTION  DOWN  TO  STATEMENT  122  COMPARES  THE  ESTIMATES 
C      OF  THE  MINIMUM  ON  THE  ERROR  CURVE 

131  ERR(KI )=ERROR(M) 
ROLDKKI  )=R1 

M  =  N 

ERRORd  )  =  ERR0R(4) 

DLTA=DLTA*VAR 

I F ( KI-1 )133,132.133 

132  KI=KI+N 
GO  TO  62 

133  TEST1=ERR(KI )-ERR(KI-l ) 
IF(TESTl) 134,135,135 

135  TEST=TEST1 
GO  TO  141 

134  TEST=-TEST1 

141  IF(TEST/ERR(KI ) -CR I T ) 71 ,71 , 136 

136  TEST=(  ROLDKKI  )  -R0LD1  (  K  1-1  )  ) /ROLDKKI  ) 
IF(TEST )137,138,138 

137  TEST=-TEST 

138  IF(TEST-VAL)71,71,139 

139  IF(TESTl) 140,133*62 

140  ERR(KI-l) =ERR(KI ) 
R0LDKKI-1)=R0LDKKI  ) 
GO    TO    62 

C      THIS  SECTION  INCREMENTS  Rl 

122  IF(L)64,62,64 

62  L=N 
RAD(M)=R1 
R0LD=R1 

63  R1=R1+DLTA 
M=M  +  N 
RAD(M)=R1 
GO  TO  79 

64  IF(M-3)65,602,602 

6  5  IF ( ERROR (M) -ERROR (M-l) ) 63 » 67, 67 
67  R1=R0LD-DLTA 

IM  =  M 

E1=ERR0R( I  M-l) 
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E2=ERRCR( IM) 

ERRCR( IM) =  E1 

ERRCR( IM-1)=E2 

RR1=RAD( IM-1) 

RR2»RAD( IM) 

RAD( IM)=RR1 

RAD( IM-1)=RR2 

M=M+N 

RAD(M)=R1 
C      THIS  LCCP  ADJUSTS  THE  VALUE  CF  ALPHA  SC  IT  WILL  PASS 
C      THE  CHECK  AT  THE  BEGINNING  CF  THE  ALPHA t  BETA  ITERATION 

79  DC  81  I=1»J1 

81  AlA(I)=AA(I)+l.E-4 
GO  TC  13 
602  IP=M-1 
C      THIS  LCCP  CALCULATES  THE  DELTA  ERROR 
DC  68  K=1,IP 

68  DERR(K)=ERRCR(K+1)-ERR0R(K) 

C      THIS  LCCP  ZEROS  THE  MATRIX  TAU»J) 
DC  69  I=lt5 
DC  69  J=lt6 

69  TA{ I tJ)=G 

C      THIS  LCCP  PUTS  ONES  IN  TA(KtK)  AND  TA(K»6) 
DC  70  K=l»5 
TA(K»6)=N 

70  TA(K»K)=N 

C      THIS  LCCP  INSERTS  THE  MARTIX  ELEMENTS  TA(K.6)  AND  PUTS 
C      ONES  IN  TA(1,K) 

DC  73  K=1»IP 

TA(K»6)=DERR(K) 

TA(K»1)=N 
C      THIS  LCCP  CALCULATES  THE  AVERAGE  RADIUS  USED  TC 
C      DETERMINE  THE  DELTA  ERRCR 

DC  73  J=2,IP 
73  TA(K,J)=( (RAD(K)+RAD(K+l>)/2.)**( J-l ) 
C      THIS  STATEMENT  SOLVES  THE  MATRIX  USING  THE  CRAM  SUBROUTINE 

B(l  )=CRAM(5# ) 
C      THIS  STATEMENT  ESTIMATES  THE  MINIMUM  OF  THE  ERRCR  CURVE 

Rl=-B(l )/B(2) 

M=M+N  I 

GO  TO  79 
C      THE  FOLLOWING  STATEMENTS  PUNCH  THE  FINAL  OUTPUT 
C      WHEN  A  CALCULATION  HAS  BEEN  COMPLETED 

71  PUNCH7 

DO  80  I =1 » Jl 

80  PUNCH1»A( I ) »B( I )»C( I )tAA(I ) 
PUNCH  8 

PUNCHl»ERR(KI-l)tERR(KI)»RltRCLDl(KI-l) 
PUNCH  101 

DC  76  I=1»J1 
76  PUNCH1»BG( I) 
PUNCH93 
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DC  500  I=1.NRCD 

FLUXR=G 

DC  501  J=1»J1 

501  FLUXR=A( J)»AIZ( Jtl )+FLUXR 
500  PUNCH99»FLUXR.I »RA(  I  ) 

DC  502  I=L1»N2 
FLUXG=G 
DC  503  J=ltJl 
50  3  FLUXG=FLUXG+B( J)*AJZ(JtI >+C ( J ) *AYZ < J» I ) 

502  PUNCH99»FLUXG»I »RA( I  ) 

C      THIS  STATEMENT  RETURNS  TC  THE  BEGINNING  CF  THE 

C      PRCGRAM  TC  READ  NEW  DATA  TC  BEGIN  ANCTHER  CALCULATICN 

GC  TC  15 
C      THE  FCLLCWING  ARE  THE  FCRMAT  STATEMENTS  USED  TC  READ 
C      IN  THE  DATA  AND  PUNCH  THE  RESULTS 
1  FCRMAT <E 16. 7 »E16.7»E16.7»E 16.7) 

5  FCRMAT(El0,4»E10.4»E10.4»El0.4tEl0.4»E10.4) 

6  FCRMAT ( I5»I5»I5»I5»I5»I5tI5.I5) 

7  FCRMAT(7X3H  A  »13X3H  B  tl3X3H  C  »12X6H  ALPHA) 

8  FCRMAT(6X6H  ERCLD10X6H  ERRCR12X3H  R112X6H  RCLD  ) 
994  FCRMAT(6X6H  RADR  »10X6H  PHIR  »10X4H  WW  ) 

993  FCRMAT(6X6H  RADG  tlOX6H  PHIG  .10X4H  WW  ) 
93  FCRMAT(5X5H  FLUX8X2H  I6X2H  R) 

95  FCRMAT(   9H  SLCP  =  0) 

96  FCRMAt(  10H  SLCPE  =  0) 

99  FCRMAT(E14.7t I6»E14.7tI6) 

100  FCRMAT<E16.7»8H  *  ERROR »E16.7» 5H  *  Rl) 

101  FCRMAT(  7X5H  BETA) 

117  FCRMAT(8X4H  SUM10X6H  SLCPE10X6H  BETA  10X6H  ALPHA) 
END 
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APPENDIX  G 


Typical  Output  Results  from 
IBM  1620  Effective  Rod  Radius  Program 


RESULTS  FOR  CASE  1C 
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1.0000000E-06 
1.1267290E-01 
4.9519509E-01 
1.0000000E-07 
9.4700000E-01 
1.1900000E-02 
41  2 
RADR 
.0000000E-99 
6.2500000E-02 
1.2500000E-01 

1.8750000E-01 
2.5000000E-01 
3.1250OOOE-O1 
3.7500000E-01 
4.3750000E-01 

RADG 
5.0000000E-01 
5.6250000E-01 
6.2500000E-01 
6.8750000E-01 
7.5000000E-01 
8.1250000E-01 
8.7500000E-01 
9.3750000E-01 
l.OOOOOOOE-OO 
1.0625000E-00 
1.1250000E-00 
1.1875000E-00 
1.2500000E-00 
1.3125000E-00 
1.3750000E-00 
1.4375000E-00 
1.5000000E-00 
1.5625000E-00 
1.6875000E-00 
1.8125000E-00 
1.9375000E-00 
2.0625000E-00 
2.1875000E-00 
2.3125000E-00 
2.4375000E-00 
2.5625000E-00 
2.6875000E-00 
2.8125000E-00 
2.9375000E-00 
3.0625000E-00 
3.1875000E-00 


4.4059408E-00 
3.6023495E-01 
3.0000000E+01 
1.0000000E-04 
2.2900000E-01 
5.6000000E-02 


PHI 
1.01153 
1.00588 
1.00014 
1.01124 
1.04386 
1.06365 
1.09502 
1.10620 

PHI 
1.15931 
1.18927 
1.21003 
1.23580 
1.24238 
1.29351 
1.29827 
1.33215 
1.36711 
1.35689 
1.36589 
1.38403 
1.36820 
1.38599 
1.39493 
1.39922 
1.39588 
1.40070 
1.43750 
1.45624 
1.47629 
1.49410 
1.47686 
1.48943 
1.50807 
1.51889 
1.47664 
1.46171 
1.42605 
1.41677 
1.40606 


R 

53E+04 

83E+04 

03E+04 

35E+04 

48E+04 

82E+04 

01E+04 

97E+04 

G 

80E+04 

87E+04 

64E+04 

13E+04 

67E+04 

79E+04 

34E+04 

34E+04 

58E+04 

59E+04 

83E+04 

01E+04 

03E+04 

86E+04 

77E+04 

22E+04 

88E+04 

22E+04 

44E+04 

05E+04 

74E+04 

96E+04 

12E+04 

00E+04 

86E+04 

22E+04 

OOE+04 

84E+04 

47E+04 

74E+04 

83E+04 


6.0638201E-03 
1.0000000E-03 
5.0000000E-03 
8.0000000E-01 
8.6000000E-04 


WW 
9.5816865E-05 
9.6392017E-05 
9.6983427E-05 
9.5973876E-:05 
9.3049491E-05 
9.1377633E-05 
8.8827885E-05 
8.7978286E-05 

WW 
8.4021765E-05 
8.1958594E-05 
8.0602182E-05 
7.8972735E-05 
7.8593324E-05 
7.5554831E-05 
7.5314384E-05 
7.3448316E-05 
7.1617831E-05 
7.2185164E-05 
7.1745558E-05 
7.0842090E-05 
7.1690302E-05 
7.0805879E-05 
7.0384294E-05 
7.0202071E-05 
7.0396815E-05 
7.0188412E-05 
6.8436376E-05 
6.7590054E-05 
6.6708813E-05 
6.5949209E-05 
6.6740481E-05 
6.6211313E-05 
6.5428052E-05 
6.4991958E-05 
6.6867407E-05 
6.7575087E-05 
6.9279635E-05 
6.9757850E-05 
7.0316495E-05 


3.2000000E+01 
5.1623768E-00 
1.0000000E-05 
1.1760000E-00 
3.8500000E-01 
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3.3125000E-00 
3.4375000E-00 

SUM 
1.2000000E-06 
3.9000000E-06 
4.0164002E+02 
1 

SUM 
6.0000000E-07 
4.4000000E-06 
4.0333814E+02 
2 

SUM 
-6.0000000E-07 
•3.8000000E-06 
4.0285732E+02 
3 

SUM 
-5.0000000E-07 
2.7000000E-06 
4.0162392E+02 
4 

SUM 
1.1000000E-06 
4.1000000E-06 
4.0266891E+02 
2 

SUM 
■1.0000000E-06 
•3.7000000E-06 
4.0255438E+02 
3 

SUM 
7.0000000E-07 
•1.0000000E-07 
4.0162228E+02 
4 

SUM 
•1.3000000E-06 
3.8000000E-06 
4.0225432E+02 

2 

SUM 
2.0000000E-07 
•2.5000000E-06 
4.0225559E+02 
3 


1.4086218E+04 
1.3980354E+04 

SLOPE 
■6.2413793E+02 
■5.2775842E+02 


7.0221602E-05 
7.0780810E-05 

BETA 
2.0958535E-00 
2.0889921E-00 


ERROR 


4.9519509E-01  =  Rl 


SLOPE  BETA 

-6.1272727E+02    2 .0958508E-00 
-5.2321984E+02    2 •0889849E-00 
ERROR    5.0019509E-01  =  Rl 

SLOPE  BETA 

-6.2379839E+02    2 •0958562E-00 
-5.3262975E+02    2 .0889994E-00 
ERROR    4.9019509E-01  =  Rl 

SLOPE  BETA 

•6.1623931E+02    2 .0958537E-00 
-5.2889298E+02    2  •0889927E-00. 
ERROR    4.9478277E-01  -  Rl 

SLOPE  BETA 

-6.2307692E+02    2.0958516E-00 
■5.2519726E+02    2.0889870E-00 
ERROR    4.9878277E-01  =  Rl 

SLOPE  BETA 

■6.2303714E+02    2.0958559E-00 
■5.3226600E+02    2.0889985E-00 
ERROR    4.9078277E-01  ■  Rl 

SLOPE  BETA 

•6.2113821E+02    2.0958538E-00 
■5.2945212E+02    2  .0889929E-00 
ERROR    4.9466681E-01  =  Rl 

SLOPE  BETA 

•6.1687242E+02    2 .0958521E-00 
■5.2661591E+02    2 .0889883E-00 
ERROR    4.9786681E-01  =  Rl 

SLOPE  BETA 

•6.2220039E+02    2.0958555E-00 
■5.3131034E+02    2.0889975E-00 
ERROR    4.9146681E-01  »  Rl 


ALPHA 
8.5293118E-02 
1.8970686E-01 


ALPHA 
8.5359468E-02 
1.8978640E-01 


ALPHA 
8.5226589E-02 
1.8962727E-01 


ALPHA 
8.5287640E-02 
1.8970030E-01 


ALPHA 
8.5340744E-02 
1.8976394E-01 


ALPHA 
8.5234420E-02 
1.8963663E-01 


ALPHA 
8.5286098E-02 
1.8969846E-01 


ALPHA 
8.5328595E-02 
1.8974937E-01 


ALPHA 
8.5243527E-02 
1.8964752E-01 
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RESULTS  FOR  CASE  1C 


SUM  SLOPE  BETA  ALPHA 

-7.0000000E-07  -6. 1875000E+02    2.0958538E-00    8. 5286 120E-02 

2.8000000E-06  -5. 2985978E+02    2 .0889929E-00    1.8969848E-01 

4.0162284E+02  =  ERROR    4.9466840E-0 1  =  Rl 
4 

ABC  ALPHA 

-5.1614880E+02  -9.912 1058E+02   -1 . 5629428E+02    8. 5286 120E-02 

1.0033036E+04  1 .77891 12E+04    3.0815911E+03    1.8969848E-01 
EROLD  ERROR  Rl  ROLD 

4.0162228E+02  4.0162284E+02    4.9466840E-01    4.946668 1E-01 
BETA 

2.0958538E-00 

2.0889929E-00 

FLUX  I        R 

9.5168880E+03  1   .0000000E-99 

9.5574740E+03  2  6.2500000E-02 

9.6797510E+03  3  1.2500000E-01 

9.8852860E+03  4  1.8750000E-01 

1.0176711E+04  5  2. 5000000E-01 

1.0557766E+04  6  3. 1250000E-01 

1.1033352E+04  7  3. 7500000E-01 

1.1609594E+04  8  4.3750000E-01 

1.2250533E+04  9  5.0000000E-01 

1.2462748E+04  10  5.6250000E-01 

1.2650846E+04  11  6.2500000E-01 

1.2819234E+04  12  6.8750000E-01 

1.2971160E+04  13  7. 5000000E-01 

1.3109092E+04  14  8 . 1 250000E-01 

1.3234944E+04  15  8.7500000E-01 

1.3350235E+04  16  9. 3750000E-01 

1.3456185E+04  17  1. OOOOOOOE-00 

1.3553797E+04  18  1. 0625000E-00 

1.3643895E+04  19  1. 1250000E-00 

1.3727174E+04  20  1. 1875000E-00 

1.3804218E+04  21  1.2500000E-00 

1.3875524E+04  22  1. 3125000E-00 

1.3941526E+04  23  1. 3750000E-00 

1.4002591E+04  24  1.4375000E-00 

1.4059047E+04  25  1. 5000000E-00 

1.4111176E+04  26  1. 5625000E-00 

1.4203434E+04  27  1.6875000E-00 

1.4281052E+04  28  1.8125000E-00 

1.4345384E+04  29  1.9375000E-00 

1.4397531E+04  30  2.0625000E-00 

1.4438401E+04  31  2. 1875000E-00 

1.4468752E+04  32  2.3125000E-00 
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FLUX  I        R 

1.4489227E+04  33  2. 4375000E-00 

1.4500378E+04  34  2.5625000E-00 

1.4502678E+04  35  2#6875000E-00 

1.4496543E+04  36  2 .8125000E-00 

1.4482336E+04  37  2.9375000E-00 

1.4460383E+04  38  3.0625000E-00 

1.4430972E+04  39  3* 1875000E-00 

1.4394367E+04  40  3.3125000E-00 

1.4350804E+04  41  3.4375000E-00 
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ABSTRACT 

The  usual  method  of  treating  the  problem  of  self- 
absorption  in  an  absorbing  media  is  by  use  of  an  effective 
cross  section  which,  when  multiplied  by  the  average  neutron 
flux  within  the  absorber,  yields  the  proper  magnitude  of  the 
flux  depression.   In  this  work,  the  problem  of  self-absorption 
was  investigated  from  the  stand  point  of  an  effective  rod 
radius.  This  approach  was  developed  by  use  of  the  P^ 
transport  model. 

The  mathematical  theory  of  the  P,  transport  model  in 
cylindrical  geometry  for  a  two  region,  absorbing  and  moder- 
ating, system  was  developed,  and  analytical  expressions  for 
the  neutron  flux  in  each  region  were  given.  Two  different 
source  conditions  were  considered  in  the  development,  and  the 
coefficients  in  the  flux  equations  were  presented. 

The  flux  depression  in  the  rod  as  the  result  of  self- 
absorption  was  approached  from  the  stand  point  of  an  effective 
rod  radius  instead  of  an  effective  cross  section  as  is  the 
usual  case.   Several  sets  of  data  for  a  one  inch  and  a  two 
inch  diameter  steel  rod  were  analyzed  by  an  IBM  1620  computer 
program  that  varied  the  radius  of  the  absorbing  rod  in  order 
to  obtain  the  best  neutron  absorption  as  calculated  using 
the  P.  transport  model. 

The  data  were  analyzed  by  fitting  the  correct  expression 
for  the  flux  in  each  region  while  requiring  that  the  weighted 
sum  of  the  square  of  the  residuals  should  be  a  minimum.  All 


of  the  data  were  analyzed  using  one  harmonic  and  the  value 
obtained  from  the  asymptotic  theory  expression  for  X    »   the 
inverse  diffusion  length,  in  the  rod.   Selected  aets  of  data 
were  analyzed  using  two  and  three  harmonics.   A  comparison 
was  made  between  P-^  transport  theory  and  diffusion  theory  by 
analyzing  the  data  using  one  harmonic  and  the  diffusion 
theory  value  for  X  in  the  rod. 

The  results  obtained  for  the  effective  rod  radius  indi- 
cated that  this  approach  to  the  problem  of  the  self-absorption 
effect  was  feasible.   However,  the  amount  of  data  used  in  this 
work  were  insufficient  to  fully  verify  the  approach,  and  con- 
siderably more  experimental  data  would  be  required  to  ade- 
quately support  the  concept  of  an  effective  rod  radius. 


